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SYNOPSIS 


OPTILLYD AITD SUBOPT BLiL COIITEOI OP IfULTL'PffiA 
LOiil) PBEQUEUCY COITTBOL SYSTEMS 
A Thesis Suhmitted 

In Partial PulfilEient of the Requirements 
Por the Degree of 
DOCTOR OP PHILOSOPHY 
by 

TEKBIiRrJA RAMilCHAHDEil MOORTHI 
to the 

Department of Electrical Engineering 
Indian Institute of Teclinology Kanpur 
September 1972 


Some computational methods to determine suboptimal 
controllers for multiarea Load Prequency Control Systems 
are presented and their nearness to the optimal controller 
is evaluated. Purther some of the realistic factors which 
are usually assumed to be absent, for the sake of simplicity, 
are considered hcze in their complete shape and their effect 
on the response- of the Load Prequency Conti-ol (LPC) system 
studied , 

The small perturbation (linear) model of a multi- 
area LPC system is considered in Chapters 2 and 3 
wherein two methods of reduction in computation are given. 

In the first method, a three-equal-area system comprising 
of 11 state variables is approximated as a tv/o -unequal -are a 
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system having 7 state variables, This is achieved by 
combining two of the areas (in the three-area system) into 
a single large area and treating it as connected to the 
third-area by an equivalent tieline. The optimal response 
of the reduced system is then determined and compared v/ith 
that of the original system. The method is then extended 
to a three -unequal-area system which is approximated as a 
tv/o -unequal -are a system. 

A second method of reduction in computation of 
optimal conti'ol of linear systems invoices the theory of 
aggregation of Aoki. In this method the original tv;o-area- 
system comprising of 9 state variables is reduced to a 
system consisting of 6 state variables which has, as its 
eigenvalues, tiie dominant eigenvalues of the original 
system. The optimal feedback controller is then computed 
for the reduced system and is used as a suboptimal controller 
for the original system. The nearness of this suboptimal 
control to the optimal one is examined by comparing the 
following, a) The suboptimal responses v/ith optimal ones, 
b) the suboptimal performance index u/ith that of the optimal 
one and c) the stability of the closed-loop systems in both 
the cases. This method involves correct knowledge of the 
eigen-vectors of tlie original system. 

In a multiarea LkC system, the tieline power flow 
between any tv/o areas is a function of the sine of the 
angular difference between the voltage vectors at either 
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end of the tie line . Por sraall perturbations, tliis angular 
difference is small and the ticline power flow can be taken 
as a function of the angular difference itself which makes 
the system linear. However, for large perturbations, the 
above assumption cannot be made and the sj^tem must be 
treated as nonlinear. In Chapters 4 and 5, advantage is 
taken of the trigonometrical sine term occuring as the 
nonlinear teim in IHC systems and the method of reduction 
by aggregation is extended to reduce a higiier order nonlinear 
system into a lower order one. The method of lulces, given 
for regulation of autonomous nonlinear dynamical systems 
is then conveniently applied for the reduced nonlinear 
system. As in the case of linear systems, ti'm optimal feed- 
back controller of the reduced nonlinear system is utilized 
as a suboptimal controller for the original nonlinear system. 
With the help of a 4th order example, it is aho\m that the 
suboptimal response closely resembles the optimal one* 

One of tiic usual assumptions that is made while 
computing optimal control of LHC systems is that there is 
no interaction between the Megavar-voltago and Megawatt- 
frequency control loops. This is not strictly true in the 
actual state of affairs. The effect of inclusion of 
excitation control on the frequency variation, which is 
obtained as a function of time, is studied in the last 
part of this thesis, A second realistic factor included in 
LHC systems and studied in this part is the construction of 
an observer in the event of nonmeasurability of some of the 
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state variables and its effect on the response of the state 
variables. It is shown that the error in the response thus 
obtained is small and also decays exponentially v/ith time, 
thus achieving a near-optimal controller v/ith the help of 
the observer, ‘The above tv/o aspects are covered in 
Chapters 6 and 7 of the thesis . 

In the author's opinion, the following contributions 
are made in this thesis; 

i) Elegant and efficient computational teclmiques are 
developed for the optimal and suboptimal control of linear 
multiarca lEC systems . 

ii) Sub optimal regulation of a nonlinear multi area 
lEC system is aichieved by reducing the order of the non- 
linear system to a computationally convenient size. 

iii) The effect of excitation conti-ol on the response 
of the lEO system is determined taking into account, the 
interaction of the area voltage on the area load also. 

iv) An observer is constructed when some of the 
state variables are assumed to be nonmcasurable and 
near-optimal control achieved by using the estimates of 

the state variables so obtained, in the optimal controller 
structure . 

Finally, concluding remarks are made on all the 
computational methods presented in the thesis. 



CHAPTER 1 


lUTRODUCTION 


The study of Load Frequency Control (LFG) of 
multiarea systems is necessitated by the importance of 
maintenance of constant frequency and tieline power 
flows in large power networks. A disturbance in any 
part of the network has its effect on the frequency of 
the entire network. This problem has been studied by 
considering a ’’coherent” group of generators tka't swing 
together under distuiliance , into a single control area. 
It 1? necessary to consider, as many areas as the number 
of such "coherent” groups. These areas are interconnec- 
ted together by means of several tielines, 

i p Q 

An LFC system has been treated in the past ’ 
as a conventionally controlled closed loop System with 
speed governor, turbine and equivalent generator repre- 
sented by their transfer functions. The system is 
controlled by a signal which depends upon the frequency 
deviation, tie line power deviation and/or their 
inte grals , 

The advent of optimal control techniques has had 
considerable impact on the analysis of LFC systems also. 
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as is tlie case with many other control systems ♦ By these 
techniques ; the closed loop controller is computed in an 
optimum fashion as a linear combination of all the avail- 
able state variables viz,, frequency deviation and its 
integral, tieline power deviation and its integral, 
governor valve position and the turbine power output^. 

In this thesis, these techniques are extended to devise 
some elegant and efficient computational methods for the 
determination of suboptimal controllers for BBC systems. 

1 .2 ERIEB ElVIBW OB PEEVIOUS WORK BOHE ON LOAD 
BHEQIENCY CONIROl 

The small perturbation model of a two-area BBC 

system was first studied on the analogue computer by 

Concordia and Kirchmayer ’ . The transfer functions of 

the governor, steam turbine and equivalent generator were 

simulated on the computer, and by taking a step load 

disturbance in one of the areas, the frequency deviation 

and tie line deviation were obtained as functions of 

4 22 

time. Bater Tan Ness ^ studied the model of a large 
system consisting of eight hydroelectric generating 
stations, each having integral control. He also presen- 
ted a method for finding the eigenvalues and their sen- 
sitivities to the parameters of the system. Hov/ever, the 
above work considered the BBC system to be linear. Bor 
the first time Aggarwal and Bergseth considered a non- 
linear multiarea BBC system with conventional controller 



3 


under condition of large perturbation of load and applied 
nonlinear Programming techniques to optimize the perfor- 
mance of the system. 

3 

Elgerd and Posha seem to have applied optimal 
control techniques for the first time to the linearized 
model of a two-area IPC system. The state variables in 
their study consisted of the frequency deviation and its 
integral, tielin^deviation and its integral, governor 
valve position and the turbine power output. They assum- 
ed that all the state variables are accessible for meas- 
urement, there is no noise in the system, and that there 
is no interaction between Hega'’'ar voltage and Megavratt 
frequency control loops. 

It is the aim of the present thesis not only to 

seek methods of reduction in the computation of optimum 

controllers for multiarea IPO systems but also to study 

the IPO system v/ithout the simplifying assumptions made 

3 

by Elgerd and Posha , 

1 .3 OUTEIIIB OP VAEIOUS CHAPTERS 

In Ohapter 2 a method is presented for the 
approximation of a three— equal-area system to a two- 
unequal-area system. The responses of the state 
variables of interest are determined in both cases and 
compared. Hext, the approximation of a three -nne qua 1- 
area system is taken up and studied in a similar manner. 
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In Ohapter 3 a method for reduction in computation of 

multiarea IiPC systems is presented. In it, the theoiy 

6 

of aggregation due to Aoki is invoked and tte original- 
system is reduced to one, of lesser dimensionality. In 
the process of reduction, it is ensured that the reduced 
system matrix has, as its eigenvalues, the dominant 
eigenvalues of the original system matrix. The optimum 
controller is then computed for the reduced system and 
is used as a sub optimal controller for the original 
system. The performance index and closed loop stability 
with suboptimal control are compared with those obtained 
by using the optimum controller. 

One of the contributions made in this thesis is 

the suboptimal regulation of the large perturbation 

model of a two-area LK! system. This aspect is dealt 

with in Chapters 4 and 5. A multiarea IPC system under 

large load perturbations is essentially a nonlinear 

system in which the tieline power constitutes the non- 

7 

linearity. Liikes has presented a method for the regu- 
lation of autonomous nonlinear dynamical systems; but it 
is not convenient to apply the same to systems of, say 

7th order and above. Here, by extending the theory of 

6 

aggregation due to Aoki to nonlinear IPG systems a]so, 
the 8th order nonlinear two-area IPC system is reduced 
to a 5th order nonlinear model. A suboptimal closed 
loop controller is then designed for the original 



system by applying the method due to Lukes to the reduced 
system. 

The effect of excitation control as well as the 
effect of load variation with voltage, on the dynamics- 
of an LKJ system are studied in Chapter 6. In Chapter 7 

O 

an observer of the Luenberger type is constructed assu- 
ming that some of the state variables are nonmeasurable 

and the responses so obtained are compared with those 

5 

obtained by Elgerd and Eoshat 

Einally, in Chapter 8, concluding remarks are 
given on the different computational methods presented 
in the thesis j also, suggestions are given for further 
research in this inexhaustible area. 



CHAPTER 2 


EEDUCTIOH OE liOiP EREaHEHCy GOHTROI SYSDEBIS 
BY APPEOXBIATIOH 


The present day power systems are quite large in 
extent and for the purpose of load frequency control, 
they have to be grouped into a number of control areas 
depending on their electrical characteristics. Eor the 
application of optimal control techniques, it is neces- 
sary to identify the state variables in each area. The 
number of such variables may be four or more depending 
upon the model chosen. Therefore the total number of 
state variables for a multiarea system may be quite 
large; and computation of optimal control for such a 
system not only presents storage problem on the computer 
but also other computational difficulties. It is there- 
fore imperative that either some kind of approximation 
is adopted or some kind of suboptimal control is computed 
to reduce the amount and time of computation. One such 
approximation technique is presented in this chapter. In 
Sections 2,2 to 2.4 reduction of a three -equal- are a 
system to a two— unequal -are a system using this technique, 
is considered; and in Sections 2.5 to 2,7 reduction of a 
three-une qual-area system is dealt with. 
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2.2 THSEE-EQUAI-AEEA STSTEK 

A three -equal -are a system is so named because all 
the three areas are of the same Megawatt capacity and 
have identical gains and time constants in the correspon- 
ding control blocks. The choice of a three -equal-are a 
system is first made for the sake of simplicity in com- 
putation, However, in Section 2.5 a more general case 
viz., a three -unequal-are a system is taken up for study. 
Eigure 2,1 gives a schematic diagram of a three-equal- 
area system. A step load disturbance is considered in 
area 1 to study the dynamic characteristics of the 
system. The three-areas are interconnected by means of 
tielines Tle^ 2 » ^i ®23 respectively. The tie- 

line angles are taken as 45® » 30® and 15° respectively, 
Hov’/ever, in a practical system these angles are deter- 
mined not only by the angles of the voltages at either 
end of the tieline but also by the angles of the tieline 
impedances. Here, for simplicity, the tielines are 
taken to be pirrely reactive and thus the impedance angle 
is exactly 90°. 


Each of the areas are represented by the state 
variables Af, AXg^, APg and representing the 

change in frequency, change in valve position as a 
result of governor action, change in steam turbine 
output and change in tieline power (in p.u.) flowing 
out of the area, respectively. In this three-area 
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system, it is enough to consider and 

because the tieline powers expressed in Megawatts, flowing 

out of the three -areas sum up to zero or 


AP 


tie 3 


= -( AP 


tiel ^rl 


+ A P. 


tie2 ^r2^ 


( 2 . 1 ) 


where P^-j , P ^2 base po;‘vers in the three 

areas. Heare, as the Megawatt capacity of the three-areas 
are equal, we have P^-j = P ^2 “ ^r3* Hence (2,1) reduces 
to 


A P , . _ = — AP, . ^ + A P , . ri ) 
tie3 1^16' tiel-^ 


( 2 . 2 ) 


It- is assumed that the load disturbance in Area 1 
is such that it results in oiiLy small variations in tie- 
line angles. Hence it is justified to take a linear 
system. Ihe block diagram of a sin^e area in a three- 
area system is given in Figure 2,2 and its mathematical 
description is given below: 

2% a 

at ‘h + h + ‘hiei = ‘ V ^ai 


at -/ifj, at) (2.4) 


P=1 

pA 


^ AP . = - ffr 

gx -gx 


1 


dt “"gi - ^ tt: ^^gi + r^_. —gvi 


ti 


(2.5) 


dt 


A Z„_, = 


1 


gvi 


I 


AX 


■gvi 


gvi ' R. ^^ci 


gvi X 


gvx 


. . (2.6) 


. .. i = 1,2,3. 
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Parameters of the ghree-area System 

The description of the parameters used in the above 
differential equations and also their numerical values are 
given below! 


t* is the nominal frequency. 

^t1 »^t2’^t3 Steam turbine time constants in Areas 1 ,2 
and 3 • 

^gvl ^^ gv2 ’'^gv3 governor time constants. 

, Eg? ^3 Self regulation coiastants of generators. 

D.J , D 2 , load frequency constants. 

^r1 = ^r2 = ^r3 = ^000 M 


‘tie 1 


’tiel 


^tie 12 ^tiel3 


^1 P 2 I 

Y ' p — Sin(5-, - ^o) X P 
•^12 ^ '^T3 rT 




Sin(*.j - 




^^tie1 ~ ^^tiel2 ■^^^tie13 


^ll 1^2 * 

= y^ --^ ~ Oos(6 ? - 6o)(A6 ^ - A 69) 

2 T ‘ ^ 


rii rj 

-X 


^*2-Cos(fi-j - «^)( A5^ -A 6^) 


13 


= T^2^/ + Ti 3(/ Af^dt -/ Af^dt) 


where 


T 


12“ X^2 Pr1 " *2^ 

= 2itr; 0.1 X Gos(5i - ^ 2 ^ 
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ana 


2ir V 


I 


13 ” 


^1i 


V. 


^T3 ^r1 


los(' 



~ 2s X 0,1 X Cos(6-j - 6^) 

Also ^ = 2ir / Af^at ana A 62 = 2s f Af 2 at 

*^tie 2 " ^^tie 3 definea in a manner similar to 
AP^i^^ ; ^23 definea in the same manner as P-] 2 * 

- fig = 45 °, fi^ - fij = 30% fig - 65 = 15° 


f* = 60 Hz 

H-^ “ Hp “ ^^5 ™ ^ sec, 
D 


-3 


= D 2 = ^ = 8*53 X 10 p.u. m/Hz 
^t1“ ^t2“ ^t3" 


T 


gvl 


= P 


gv2 


= 0.08 sec. 


R 


.-j = Eg = = 2 ,4 Hz/p .u. m 


Also 

. 2H. 

■^kpi ^ ^ ®pi f* ’ i = 1>2,3 

The system descrihea in (2.3) to (2,6) can be expres- 
sea aynamically as the vector aiffexential equation 


x = Ax + Bu+rA^ (2.7) 

Here x is an 11 -element state vector consisting of the 
state variables Af-j , Afg, Af^j A , AX^g, AX^^, 

^ ^g1 ’ *^g2’ ^^g3' ^^tiel ^^tie2’ ^ ^ 

aisturbance vector consisting of AP^-j , ^^^2 ^^d5 


its elements 
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Por applying the optimal control theory, it is 
necessaiy to eliminate the third term in the right hand 
side of (2.7) thus reducing the equation to the familiar 
form: ‘ 

X = A X + B u (2.8) 

This has been achieved by redefining the states and 

controls in terms of their steady state values occurring 

3 

after the disturbance is applied . 

The performance index to be minimized is taken as 

J = i / (x^ Q X + u^ E u) dt ■ (2,9) 

Here R is taken to be a 3x3 identity matrix. The matrix 
is selected such that the frequency deviation from the 
nominal value is minimum in the three areas; also the 
tieline power deviation from the nominal value should be 


and ^ are described 

T 

0 0 0 0 

0 0 0 0 

0 0 0 0 


minimum for the three tie lines. 


The system matrices B, r , A 


below: 


B = 


0 0 

0 0 

0 0 


0 

0 

0 


T 

gvi 

0 

0 


0 


1 


T 


gv2 
0 


0 

0 


0 


0 


T 


gv3 


( 2 . 10 ) 



u 


r = 


A=: 


~1 

"Pt 

0 


0 

-1 


0 

0 

0 

0 

0 


P 0 

h 









T 
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0 

0 0 

0 

0 

0 0 

0 


t* 

2H2 

0 

0 

0 0 

0 

0 

0 0 

0 


0 


f* 

2h:x 

0 

0 0 

0 

0 

0 0 

q_ 




0 





* • 

(2.11) 

0 

0 

0 

0 

0 

^ m 

0 

0 

on 






•^pt 



'^pT 

'' 

-1 

0 

0 

0 

0 

0 

■Scp2 

0 

0 

'■\p2 

^p2 

p2 

-®p2' 








0 

-1 

‘^-■Z 

P3 

0 

0 

0 

0 

0 

4 

tm. 

P3 

2p3- 

®P3' 

0 

0 

-1 

0 

0 

0 

0 

0 

0 

0 

gvi 

■1 









-1 


0 

-1 

0 

0 

It 

0 

0 

0 

0 



0 

-1 

_ 

r\ 

-1 


0 


0 

0 
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T • 
gv3 

0 

0 

0 

0 

1 

0 

0 

-1 

^tT 

0 

0 

0 

0 

0 

0 

0 

1 

^t2 

0 

0 

-1 

^t2 

0 

0 

0 

0 

0 

0 

0 

1 

'P 

-^t3 

0 

0 

-1 

"t3 

0 

0 


-T 

13 ° 

0 

0 

0 

0 

0 

0 

0 

12^25 

*^23 ^ 

0 

0 

0 

0 

0 

0 

0 


( 2 . 12 ) 



15 


q = 


1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 


... (2.13) 

Pontryagin’ s minimum principle (optimal control theory) 
is applied to find the closed loop controller of the form 
(Appendix B) 

u(t) = -K x(t) = P x(t) (2.14) 

where P is the solution of the matrix Riccati equation 

(2.15) 


P = -P A - A^ P + P B R”"' B^ P - Q 


and is a 11x11 matrix. In riew of infinite time involved 
P becomes zero and (2.15) becomes simply a matrix quadratic 
equation given by 

-1 m 

(2.16) 


0 = -P A - A^ P + P B R”^ B^ P - 'Q 


Ihe following procedure is adopted for obtaining the 
solution of this equation. At first the differential 
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equation (2.15) is solved for 8 or 10 seconds such that 

the matrix P so. obtained makes the closed loop system 

wl T 90 

matrix (A - B R B P) as a stable matrix . Then this 

result is fed as initial solution to the algorithm given 
11 1 2 

by Man ’ and the correct solution matrix P is obtained. 
Appendix C gives a brief description of this algorithm of 
Man, The closed loop' system matrix now becomes: 

X = (A - B R"*^ B^ P) X , (2.17) 

and if a load disturbance is applied in Area 1, this 
becomes 


X = (A - B R”^ B^ P) X + (2.18) 

where ^ is a 11x3 load disturbance distribution matrix 
as described in (2.11). The vector 



^d1 


~o.of 

c> 

11 

^d2 

= 

0.0 


_^d3_ 


0.0 


( 2 . 19 ) 


Equation (2.19) sho7/s that a load disturbance of 
0.01 p.u. of po7/er is considered only in area 1 . The 
initial condition vector is taken as x(o) = 0. The dynamic 
equation (2,18) is solved and the responses of the state 
variables of interest viz., ^f-j » ^^tiel and A f 2 

are obtained. As the disturbance is in Area 1 , the tie- 
line power interchange error of Area 1 , ^^tiel’ 
imp ortance . 


is of 
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2.3 APPEOXBIA.IE TWO-AEEA SYSTEM 

In Section 2.2, the three-equal -area system is 
studied in its complete shape without any approximations. 

In this section reduction in computation for the three- 
equal-area system is achiered by app 2 ?oximating the same as 
an equivalent two-area system. Here, the Areas 2 and 3 
are combined together to form a sin^e-area (Eigure 2,5). 

As we have considered the three-areas to have identical 
parameters, the capacity in Ifflf of the combined system 2 and 
3 will be double that of any one of the three areas. The 
tieline here consists of two branches Tie-j 2 and Tie^^. As 
our ultimate aim is to consider this system as a tv/o-unequal 
area system, the tielines lle -|2 ^^®13 combined 
together into a single equivalent tieline connecting the 
two-areas. The parameter T -^2 this new system is approxi 
mated as the sum of the parameters T-j 2 (old) and T-jj(old) 
corresponding to Tie-j 2 and respectively. Thus 

T^2(^ew) = T^2 (o 1^) + l'i3(old) (2.20) 

or 

T-j 2 (iiew) t= 0,1 X CoslS*^ x 2^^ + 0,1 x Cos30o x 2'ff 

.tv 

= 0.1 X 2Tfx (Cos 300 + Cos 45°) (2.2l) 

This may also be interpreted as taking the new 
tieline to consist of two component tielines, the angle 
for both of which is the same and is approximated as the 
average of the angles of the corresponding tielines of the 
t hre e -are a-sys tern. 




Area 2 


Fig. 2*3 Equivalent two-area system 



The relationship between A P, . „ and AP^. . 

tie<i tie I 


becomes 


AP^.. 


= a. 


AP^_. 


( 2 . 22 ) 


■tie2 “ nz “^tiel 

’Where ^^tTe2 unit of power. As 

their Wff equivalents are equal and opposite in sign, 


^ ^tie2 ^r2 "■ ^ ^tiel ^r1 


(2.23) 


where P^-j and P ^2 ^re the base powers in the new Areas 1 
and 2 respectively. So 


^^tie2 “ ^^tiel 


(2.24) 


Because we have combined the areas 2 and 3 of the three- 
equal-area-system, the base power of Area 2 in the new 
configuration vd.ll be twice that of either one of the 
areas . Thus 


^r2 ^^rl 


Prom (2.22), (2,24) and (2.25), we get 
P . 

^2 = -(p^) = 

^f2 


(2.25) 


(2.26) 


Because of this relation, we need not consider ^ ^ie2 as a 
a state variable, Por'the equivalent tv/o area system, the 
system vector differential equation becomes 

X = A X + B u (2.27) 

where x is a 7-element vector consisting of the state 
variables A f ^ , Afgj AXg^-[,AXg^, t ^ ^g 2 ^^tie1 * 
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The matrices A(7x7) and B(7x2) are given by 



(2.29) 


All the parameters used for the elements of the matrices A 
and B are as given in Section 2.2 except for the value of 
a-j 2 which is given in (2.26). 

■ The performance index to be minimized is 

-f 00 m m 

/ (2£'^x + uEu)dt (2,30) 

o „ 

The matrix is selected taking into consideration the 
fact that the state variables A f ^ , Af 2 (new) and ^ 
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should have minimum variation, Thus 




0 

0 

0 

0 

0 

0 


0 0 0 0 0 
1 0 0 0 0 
0 0 0 0 0 
0 0 0 0 0 
0 0 0 0 0 
0 0 0 0 0 
0 0 0 0 0 


0 

0 

0 

0 

0 

0 

1 


(2.31) 


R is taken as a 2x2 identity matrix. Pontryagin‘s minimum 
principle (Appendix B) is applied to get the control vector u 
as 

u = P X (2.32) 


where P is the solution of the matrix Riccati equation 

P = -P A - A^P + P B R””* B® P - a (2.33) 

Here P is a 7x7 matrix. The responses of the state 
variables & f-] » * ^2’ ^ ^tiel determined by solving 

the closed loop vector differential equation with a step 
load disturbance of 0,01 p,u. of power in Area 1, Por this 
purpose, (2,27) is augmented v/ith load disturbance as 

X = (A - B R"' B^ P) X + (2,34) 

Here r is a 7x2 load-disturbance-distribution matrix. 


r 



0 


f* 

2H2 


0 

0 


0 0 0 0 

0 0 0 0 


(2.35) 
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and 



^^d1 


0.01 

_^^d2_ 


0.0 


( 2 . 36 ) 


Equation (2,36) shov/s that load disturbance is considered to 
occur in the same area in both the three-area and two-area 
cases. Thus the responses obtained in both cases can be 
compared to give an idea about the efficacy of the approxi- 
mation. 


2.4 GCaiPUT^IOH'AIi EESUIffiS FOE THREE-EQUlIi-4EEA SYSIEM 


(a) Ihree-area System filth order) 

The B matrix (2.t0)is computed as 
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The responses of the state variables Af^, a f 2 are 
shown as continuous curves in Figure 2,4 and that for ^~tie1 
is shown as continuous curve in Figure 2,5* The response of 
A f^ is found to be almost identieal to that of A f 2 and hence 
not given. 


b) TwO"-area System (7th order) 


The 

A matrix of (2,28) is computed 

as 
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0.0 

0.0 

0.0 6.0 

0.0 

- 6,0 

0.0 

-0.05 

0.0 

0.0 0,0 

6.0 

3.0 

' -5.208 

0,0 

-12.5 

0.0 0.0 

0.0 

0,0 
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0 
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The 

B matrix of (2.29) is computed 

as 
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— ,T 
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The closed loop matrix (A - B B^ P) of (2,34) 
is computed as 
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time, sec- 



-Response in 3 area system (t!XI1| 

Response in 2 area system (7X7i| 


Fig. 2-4 Response of Af, and A% 
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The responses of the state variables a f-| of the 
smaller area (Area 1) and ^±2 of the larger area (Area 2) 
are given as dashed curves in figure 2,4. AlsoAP^^^^, the 
total tieline power change out of Area 1 is given as dashed 
curve in Pigure 2,5. 


2.5 0:HEEE-UlEaUAI.-ASSA SYSTEM 

As is evident from the title, here the three-areas 
are considered to have different base powers and correspon- 
dingly different parameters. However for simplicity of 
computation, per unit values of the parameters like inertia 
constant H, load frequency constant D, self regulation 
constant R are all considered to be equal in the three areas. 
There is no loss of generality in this v/ay of treatment 
because what is more important is consideration of different 
base powers in the three areas. 


In Eigure 2.6 is given a schematic diagram of a three- 
unequal-area system in which the capacity of Area 2 is 
1 .2 times the MW capacity of Area 1 and the I/IW capacity of 
Area 3 is 1 .5 times that of Area 1 . The tie line angles for 
the three tielines Tie-jg, Tie-^^ and T 1 q 2 -^ are taken as 45®, 
30° and 15° respectively as in Section 2.2. The state space 
analysis of this 11 state variable system is the same as in 
Section 2,2 except that the differential equations pertaining 
to A f j and changed as a consequence of the 

different MW capacities of the areas . 
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(i) A 


^ it ^^3 “ “^3^ ^3 ■" ^^tie3 ^^3 "*^<13 

Considering a =0.0 and rewriting the above 


equation, 

It *^3 "" *■ ^aH^“ ^^3 “■ 2^ '^^tie3 ■*■ 2^ ^^g3 


( 2 . 38 ) 


Using the fact that the total tieline povrer in Elf flowing 
out of the three areas is zero, 


‘hies h3 = -(‘hlsl hi +‘hle2 ha^ 


(2.39) 


‘hle3 = -^‘hle1 ‘Wig) 


Here 


( 2 . 40 ) 


""rl 1 


^r2 _ 1 .2 

^r3 ^ '-5 


Hence 


*hle3 “ “ hle2 T^) 

Substituting this value for ^^tie3 ^ 


(2.41) 


1.5- 


. f*P, f^ ^ X 1 

dt ^^3 " “ -2H^ ^3 2H^ ^^g3 2H^ .1.5- 

^^tie2 ^ lit (2.42) 

Phus the elements A(3, 10) and A(3, 11) of the 11x11 system 
matrix A in (2.12) are modified as 

Af5 101 - Xt 1 _L _ _ _L 

-fl-t3,io; - 2JJ^ X X 
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and 

A(3,11) 


TT ii2 _ "kp3 ^ 

2H3 1 .5 " Tp3 1 .5 


the rest of the elements of the 3rd rov/ remaining the same. 


‘hie2 = 


®23^ ‘^2 " 

Here 

H^-i -I 

^ (2.44) 


^ 2 l " ^12 ^ 5^2 “ ®12 ^ t .2 


where 11-^2 As expressed in p.u. of base pov?er • ^23 

remains the same provided it is expressed in p.u. of P^2 
only. Hence the elements A(1.1, 1) and A(1 1 , 2) of the 
matrix A in (2.12) are modified as 


and 


A(11, 1) = -.P-,2 
I 

A(11,2) = (1^2 ^ T^2 + '^23> 


(2.45) 


( 2 . 46 ) 


the rest of the elements in the 11th row remaining the 
same as given in (2.12). Also all the elements in the 
other rows except 3 and 11 remain unchanged. 


In this three-^mequal-area-system the following 
values are assumed for '^^2^ ^13 ^23' 

1^2 = X Cos 45° X 2 IT expressed in p.u, of P^-j 

1-^3 = 0,1 X Cos 30° X 2 TT expressed in p.u, of P^^ 

^23 ~ ^ ^ ^ expressed in p.u. of P^2* 

The matrices B, E and also the theory and analysis remain 
the same as in Section 2.2. 
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A step load disturbance of 0.01 p.u. of power is 
considered in Area 1 and the responses of the state vari- 
ables of interest viz. Af^, A f ^ and ^ obtained, 


2.6 APPHOXBIAIE ITfO-UlE QUAL-AREA SYSTEM 


Here t?/o of the areas in the three -are a system 
are combined together as is done in Section 2.3, thus 
converting the tliroc-area system into an equivalent tvro- 
arca system. The system matrix A is having the same ele- 
ments as in (2.28) except for the following slight modifi- 
cations. 


As we have to combine the Areas 2 and 3 as a e cmpo- 
sitearea, their combined capacity (and inertia) is 2,7 times 
the capacity (and inertia) of Area 1 . Hence 

1 


A Pi 


Map 
P^2 tiel 


■tie2 P •- ” c? .7 (2.47) 

In the differential equation pertaining to A f^ the value 
of 3^2 taken as (-1/2.7) in viev/ of the above relation. 


Thus only the element A(2,7) of the matrkc A in (2.28) is 
changed, the rest of the elements remaining unchanged. 


The matrices B, E as v/ell as the theory and 
analysis remain the same as given in Section 2,3. 

A step load disturbance of 0.01 p.u. of power is 
considered in Area 1 and the responses of Af|, Af 2 (new) 
and 


are obtained as before 



2.7 COMPUTATIONiOj RESULTS ECR THElEB*-UNBCiUAL-ARSA SYSTEM 
(a) Three -Unequal -Are a System (11th orderl 

The cloSBd loop matrix (A-EK), 11x11 in this case is computed as 
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Figure 2.7 and that obtained for A is shown as con 




Rt-sponse m 3 area sysfem (t! X llj 
Response m 2 area system 17X7) 
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(b) Equivalent Two -Are a System 

The closed loop 7x7. matrix (A-BK) for the equivalent 
two- area system is- computed as 


-0,050 

0.0 

■ 0.0 

0.0 

6 ,0 

0.0 

6 ,0 

0,0 

-0.05 

0,0 

0.0 

0.0 

6 ,0 

2.222 

-12.403 

-0.906 

-15.339 

0.186 

-11 .859 

0.462 

13-025 

1.654 

-15.050 

0.185 

-1 5 .844 

1 .274 

-14.222 

-1 .360 

0,0 

0.0 

3.333 

0.0 

-3.333 

o 

• 

o 

0.0 ■ 

0.0 

0.0 

0.0 

3.333 

0.0 

-3 .533 

0.0 

0.989 

-0.989 

0.0 

0.0 

0.0 

o 

o 

0.0 


The responses obtained for and Af2 in this case are 
shown as dashed curves in Eigure 2,7 and that obtained 
for A is shown as dashed curve in Eigure 2.8, 


2.8 COHGLUSIOITS 

In this chapter, a three -e qual-area system is first 
taken up and approximated as a two-unequal-area system. 
Optimal controllers are constructed and responses obtained 
for A f ^ ^2 ^^tie1 both cases. 'Eigures 2,4 and 2,5 
show that the responses obtained -with the two-area repre- 
sentation are very near to the optimal responses. Hence 
it can be said that the two-area approximation is quite 
good. The study of a three-une qual-area system and its 
equivalent representation is next taken up. The responses 
obtained in this case and given in Eigures 2,7 and 2,8 
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alao ooafira the abovo results * It is soon f rom tho above 
four figures that there is a saall error in the steady state 
valus in each of the responses obtained with the two-area 
reprosentetions . 

fho study saad® in this chsa>ter indicates tJmt the 
approxlaatlon can be done, for four-area and other hi^er 
order aulti-aroa systems with good results. It is necessary 
that the area in idiioh the disturbaxuse occurs , be singled out 
as a separate area in tte reduced system also, as otherwise 
the, ijppropciiaation nay not prove to bo advantageous. 

In the three-area as well as the equivalent two-area 
systems studied above | integral feedback is absent, i.e,, 
the states ooxEresponding to / &fdt, / Al^^j^^dt of the indi- 
vidual areas aro absent. Hence the feedback control in all 
the oases I is unable to bring back the response of the state 
variables to zero stea^ state* Evidently, If integral 
feedback is also present, the responses ultimately go to 
sero steady state, Henoj there will not be any steady state 
error and tho approixination should prove to be much better, 

Ihe &mriiig in computational effort as well as In 
ooaputational time accruing out of suoh approximation is 
oonsld© rabl® , Ihis could be seen from tho fact that for a 
11th order system (three-area system) matrix liccati equation 
of 66th order is to be solved; Mdiexeas for a 7th order system 
(equivalent two-area system) the order of the matrix Eiccati 
equation to be solved is only 28, 


CHiPTER 3 


SUBOPTB/IAL COM}ROL USING AGGESG^JTION 

In Chapter 2, reduction in computation of linear 
multiarea LEO systems Mas achieved by approximating them 
as lesser oz’der multiarea systems. Here a different 
approach is attempted with the same goal in mind, Eor 

g 

this purpose, the theory of aggregation due to Aoki is 
applied. 

In Sections 3.2 and 3.3 the theory of aggregation 

g 

and suboptimal control by aggregation are reviewed beiefly. 
Sections 3.4 and 3.5 deal with a method of computing optimal 
and suboptimal performance indices. In Section 3.6 deter- 
miration of closed loop stability is discussed; and finally 
the computational results are presented in Sections 3.7 to 
3.10. 

3.2 IHSORY OR AGGREGATION^ 

Let the dynamics of the system which is linear, be 
given by the vector differential equation 

x = Ax + Bu (3.1) 

Here z is an n-vector 
u is an r-vector 
A is an nxn matrix 


B is an nxr matrix 
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Let ^ (an m vector, m < n), called the aggregated 
state vector be related to x by the relation 

z(t) = Ox(t) (3 .2) 

where C is an ircai matrix: and its m rows are the m dominant 
eigenvectors of A . It is assumed that the rank of 0 = m. 

Assuming certain matrix equalities (given below) to 
be satisfied, the response of the above system can be appro- 
ximately determined knowing the optimal feedback coefficients 
for the reduced system, viz., 

z=Pz + G-u (3 .3) 

with the assumption that z(o) = C x(o). 

Here J and G- are related to A and B by 
B 0 = C A (3.4) 

and 

G = G B (3.5) 

If A and G satisfy the matriic equation 

G A = G A G^(G G^)*"'' G (3.6) 

then B in (3.4) is given by 

B = G A c'^(C G^)”"' (3.7) 

The aggregated state vector z satisfies the dynamic 
equation (3.3) vd.th B given in (3.7). 

Here G is called the aggregation matrix: and B, the 
aggregation of A. G is the important design parameter in 
constructing the reduced or aggregated system. The choice 
of G is to be made in such a v/ay that error in modelling 
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the original system by means of the reduced system is 
minimized, lo accomplish this, it is necessary that the 
matrix G be selected such that its m rows are the m eigen- 
vectors of tile transpose of the A matrix, corresponding to 
the m dominant eigenvalues. Moreover, P is chosen such that 
it retains the dominant eigenvalues of mathematically 
this means that S' should satisfy (3.7). 

3,3 APPlilCAllClT 10 SE'lSrillilS GUBOPlBiil CONTROL 

let the performance index of the original system be 
given by 

00 m rp 

J=j(x<Qx+uRu)dt (3.8) 

o - - 

wliere -Q . is positive semidefinite and R is positive defi- 
nite, The reduced mth order system is given by (3.3) v/ith 
the initial condition z(o) = C x(o). 

let u = -K z be the optimal control policy for the 
reduced system where 

K = R-’* P (3.9) 

and P is the solution of the algebraic matrix Riccati 
e quation 

P^P + PS'-PGR”'’ G®P + <^ = 0 (3.10) 

Here is defined as in (3.12) below. The smaller system 
is considered to have the following performance index with 
as the weighting matrix for the reduced states. 
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0© m ® • 

“ / (a %i £ + !i B a) 


( 3 . 11 ) 


Where the appropilatt choice of ^ la 

= (c c’')"'o a 0 0®)"’ 


With this choice of 

» " % ; ' 

% ^ If Jig 

% « ■•'•K g w -JE 0 


( 3 . 12 ) 


( 3 . 13 ) 


beooH®® the auboptlnal control, for the original system, 
whioh yery nearly approx: Imates the optimal control. With 


the abote ettboptijoal controlf the original system is gowr- 





"the equation 


(A — B K C) ^ 


(3.14) 


fhe optimal feedback system is given by 

JK* » U - B K») ^ (3.15) 

'«» *> , I |r H# #» ## #> 

In (3*15) tba optimal control is given as 

a -K* 2f (3.16) 

where 

K* » 1"* B 0?* (3,17) 

and I* satisfies 

A® f * + I* A .. T* B B® f* + Q » p (3.18) 

3 A MEASUEB Of BBEfOEMAHOB IICDEX flfH AND 

SUBOBJIMAIi CONffiOBS 

®he following is a brief summary of the results of 
levins and Athens who have given a method for obtaining 
a measure of tlie perfomaance index with optimal and sub- 
optiiiial oonti'ols# 
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Ths closed loop system, after incorporating the 
sub-optimal control gi-^ren by (3.15) becomes 

i(t) = (A - 3 K 0) x(t) (3.19) 

Ihe solution of this differential equation becomes 


x(t) = 0i'b} O) x(o) (3.20) 

vifiiere )Z^(t, O) = exp([A -.BKClt) (3.21) 

Substituting (3.20) into the performance criterion given 
in (3.8), 


J = x^(o) 


i / iZf^(t, 0) (Q + cVEE:0)^(t,0)dt 


x(o) 

( 3 . 22 ) 


Equation (3.22) shows that the performance criterion depends 
on K, C and the initial condition x(o). We can eliminate 
the dependence on x(o) by averaging the performance obtained 
for a linearly independent set of initial states which is 
equivalent to choosing the initial state x(o) to be a random 
variable uniformly distributed on the surface of the n- 
dimensional unit sphere. Then the expected value J of J 
becomes 


J = 


1 _ 

2n 


QD 

/ tr 


j2J^(t,0)(Q+cVBKC )j2J(t,0) 


dt (3.23) 


J is now independent of the initial state and can also 
be interpreted as a completely determinlatio performance 
criterion in which the control consists of product of K 
and 0 and the state is j2^(t,0) the fundamental transition 
matrix. In (3.23) we can remove n which is only a constant. 
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The optimization problein^ow becomes: 
G-iYen the performance index 




J = 


o 


/ tr 


/(t,0)(4it.oVEKC) 


dt 


(3.24) 


T'ind -IC’^hich minimizes the performance criterion 
( 3 . 24 ) subject to the constraint imposed by the system 

j^(t,0) = A ~ B K i2^(t,0) , 0 ( 0 , 0 ) = I (3.25) 


The main result of a theorem proved by Levine and 

1 5 

Athans' for such a problem is as follows; 

Let 


H = (A - B K 0) 

(3.26) 

Assuming that H is stable, then in order for K to 

be optimal 

for the original system, it is necessary that 


K = £*“'* B^ T V 0^ [p 

( 3 . 27 ) 

wheiB 


00 m T* I! 

T = / exp(irT ) (Q + O K EE:C)exp(HT )d t 

(3 .28) 

^ JI 

V = / exp (Ha ) exp (ii ) d® 

( 3 . 29 ) 


o 


Alternatively, assuming tliat K, T, and Y are solutions of 
( 5 . 27 ) to ( 3 . 29 ) respectively then T is also a positive 
definite solution of 

0 = TH + H^T + <5+C^K^EKC (3.30) 

and V is a positive definite solution of 
0 = VH^ + HV + I 


(3.31) 
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Also iz assume that exists and denoting T under 
this condition as 1* (3.2?) reduces to 

K = R"’’ T* C"'' (3.52) 

arid (3.30) reduces to 

0 = T*A + A^ 1* + Q - T* B R’’' B^ T* (3.33) 

which is the same as- (3.18). 

3.5 BiBGRilDAlIOR IE PERP0H.1ABCE Iin)EX WITH SUBOPTEvIAL OOilTROI 

Prom (3.24), (3.28) and (3.30) it is clear that the 
trace of the matrix T is a measure of the performance using 
suboptimal control given by (3.13), irrespective of initial 
conditions j also the troce of the matrix T* gives a measure 
of the performance index using optimal control. 

Por the load frequency control problem under study, 
the matrix T as given by (3.30) for the tv/o cases of aggre- 
gation considered (viz, 6th order and 4th order), and the 
matrix T* as given by (3.18) are computed. The traces of 
the matrices T and T* are then computed and compared to 
assess the degradation in performance index, 

3.6 STABILITY OP TIE 0 IDSEB LOOP SYSTEII WZDH 3U30PT BIAL 
C aWROL 

The stability of the closed-loop system using sub- 
optimal control is determined by directly cenputing the 
eigenvalues of tiie closed-loop matrix (A - B K C) for the 
two cases of aggregation. However, this could also be 
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2 1 

detennined by liapunoT's equation. If the solution matrix M, 
v/iiich is tiie solution of the matrix equation 

(A - B K 0)^ i: + M(A - B K C) = -IT (3.34) 

is a positive definite matri^c v/ith the matrix IT being ar- 
bitrarily chosen to be positive definite, then the closed 
loop system matrix (A - B K C) is stable. 

3.7- GOlIPUlAlIOrT 01 SUBOPIB-IAL COiaBailEES BOB 61H OPDER 
AGGEEGATIOIT 

Ihe A and B matrices for the load frequency control 
system are as defined in (A. 8) and (A. 9) of Appendix A and 
computed as 
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Here the Matrix G = 


0.0 

0.544 

0,213 

0.384 

0.102 

0.544 

0.213 

0.384 

0.102 

0.875 

0.0 

0.0 

0.0 

0.0 

-0.372 

-0.146 

-0.263 

-0.070 

0.0 

0.447 

-0.108 

-0.360 

-0.100 

-0.447 

0.108 

0.360 

0.100 

0.0 

0.264 

0.261 

0.106 

0.0 

-0.264 

-0.261 

-0.106 

0.0 

0.0 

0.0 

0.144 

0.588 

0.175 

0.0 

0.144 

0.588 

0.175 

0.0 

0.0 

-0.291 

-0.132 

0.0 

0.0 

-0.291 
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0.0 


T 

The six rows of C are tte six eigenvectors of the matrix A 
corresponding to its six dominant eigenvalues. Also here 
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0 
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0 
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0 
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0.0 

0,0 

0.0 
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0.0 

0 

• 

0 

The 

solution 

of tile 

matrix 

Riccati equation for the 

reduced 


system is the 6 x 6 matrix 


2.790 

1 .109 
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Ac'iortiiiiglyi K tor tliia 6th reduced model given by 

'''In' 

I * rs 

is computed I tlien the astern matrix incorporating the eub- 
optlffiX feedback becomes 

(A - B K 0) « (A - B G® P C) 

She ears® is computed as 
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-0.050 

6.0 

0.0 

3.270 

0.0 

0.0 

0.0 

o 

» 

o 

0,0 

0,0 

-3.333 

5.333 

0*0 

0.0 

0.0 

0.0 

-8.639 

-4.200 

-15.190 

-15.190 46,139 

-8.299 

-2,428 

-2 .207 

—0 *445 

0.0 

0.0 

0.0 

0,0 

0,0 

0.0 

1.0 

0.0 

0.0 

0.0 

3,270 

0,0 

0,0 

0.0 

-3.270 

-0.050 

6,0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0,0 

C.O 

-3*533 

3.333 

8.639 

-8.299 

—2.428 

-2.207 

-0,445 

-4.200 

45 *1 90 

45.190 

-16.139 


Ite responoe of tlie feedback s^tem which is governed by the 
differential equation 
s» ( A « B 1 0 ) 

is dete mined for a load disturbame of 0*01 p.u, of power 
in the fiiwt area* 

3.8 COIOTMSICM OP SUBOPS IMAL OOSSEOmSS PCE 4SH QEDBE 
AGGffiSaATICN 

Her© the matrix! 0 is computed as 
0.0 0.544 0.213 0,384 O.m 0.544 0.213 0,384 0.102" 

0.875 0.0 0.0 0.0 0,0 -0.372 -0.146 -0.263 -0.070 

0.0 0.447 “0.108 -0.360 -0.100 -0.447 0.108 0.360 0,100 

0.0 0.264 0.261 0.105 0.0 -0,264 -0.261 -0,106 0.0 
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The four ravs of 0 are the four eige sectors of the niatrioc A 
corresponding to its four dominant eigenvalues . 


0.0 0.0 0.0 0.0 

0.0 0.0 0.0 0.0 

0.0 0.0 -0.497 3.522 

0.0 0.0 -3.522 -0.497 


1 .281 1 .281 
0.0 -0.875 

- 1.248 1.248 
0.0 0.0 


The solution of the matrix Eiccati equation for the reduced 
system is the 4r4 raatrix 


0,836 

1 .1.09 

0,173 

0.145 

1 .109 

3.247 

0.506 

0.425 

0.173 

0.506 

0.980 

-0.269 

0.145 

0.425 

-0.269 

1 .499 


-1 T 

Here also K = E G 1* is computed; the closed loop system 
matrix given hy (A - B K 0 ) = (A - B E"*”^ P G) is canputed 


0.0 

0.545 

o 

• 

o 

0.0 

0.0 

-0.545 

0.0 

0.0 

o 

• 

o 

0.0 

1 .0 

O 

• 

o 

0.0 

0.0 

0.0 

0.0 

0.0 

-3.270 

-0 .050 

6 .0 

0.0 

3.270 

0.0 

0.0 

0.0 

0.0 

0.0 

-3 .333 

3.333 

0.0 

0,0 

0.0 

-8 .639 

-1 .929 

-10.546 

— 9 .500 

-14.846 

-6 .028 

2.217 

3.683 

o 

• 

o 

0.0 

0.0 

0.0 

0.0 

0.0 

1.0 

0.0 

0.0 

3.270 

0.0 

0.0 

0.0 

-3 .270 

-0.050 

6 .0 

0.0 

0.0 

0.0 

0.0 ' 

0.0 

0.0 

0.0 

-3.333 

8.639 

-6.028 

2.217 

3.683 

0.848 

- 1 .929 

-40.546 

-9.300 
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As in the case of 6tli order aggregation, the response of 
the closed loop sj^stem governed by the differential equation 

X = (A - 3 K C) X 

is determined for a load disturbance of O.Ol p.u. of power 

in the first area. The optimal closed loop response obtain- 
3 

ed previously for the original 9th order system is used for 
comparison. 

The optimal and suboptimal responses for the cases 
of 6th and 4th order aggregations for each of the state 
variables of interest are drawn on the same figure to give 
an idea of the nearness of the suboptimal response to that 
of t}» optimal one . Figures 3.1, 3.2 and 3.3 are drawn for 
the state variables viz. frequency deviation Af^ in Area 1, 
frequency deviation Af2 in Area 2 and tieline power error 
Altiel respectively. 

3.9 00LI3DTATIOT OF DEGRADATIOl'T III PEEFOffl.IAiTCE 

The trace of the matrix T* of (3.18) obtained for 
optimal control is computed as 8.8053. For the suboptimal 
control by 6th order aggregation, the trace of the matrix T 
given by (3,30) is computed as 8,9189, Similarly for ths 
4th order aggregation the trace of T is computed as 9.5608. 
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3,10 S2ABILITI OP THE CLOSED LOOP SYSTEM - CCS'/I?UT>TIOiT 
OP EIGEIiYLLUSS 

The eigenvalues (some of which are complex) of tlie 
closed loop sysiem matrix (A - B K 0) in the case of 6th 


order agg 

i'regation are computed as 

1) 

-0.505024795 


O 

« 

O 

2) 

-0.932741128 

+ 

jo.o 

5) 

-1 .623612400 

-f 

30.0 

4) 

-1 .842025460 

+ 

33.731263460 

5) 

-1 .842025460 

- 

33.731263460 

6) 

-2.872359200 

+ 

33.451757180 

7) 

-2.872359200 

- 

3 3. 451 7571 80 

8) 

-13.265108300 

+ 

0 

• 

0 

9) 

-13.290158100 

■f 

30.0 


In the case of 4th order aggregation the cl osed loop 
eigenvalues are: 


1) 

-0.505024795 


30.0 

2) 

-1 .296577550 


32.512660310 

3) 

-1.296577550 

- 

32.512660310 

4) 

-1 .497734720 


30.0 

5) 

-1 .623612400 

-h 

30.0 

6) 

-1 .842025460 

4“ 

33.731263460 

7) 

-1 .842025460 

- 

33.731263460 

8) 

-13 .265108300 

+ 

3O.O 

9) 

-13.290158100 


3O.O 
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Ste Oigonvalues of tte open loop ayoteu aatrlx A of (3.1) 


also given belov? for 

' comparisons 

1) 

0.0 

+ 30.0 

2) 

0.0 

+ 30.0 

3) 

"0 .497296223 

+ 33.522088710 

4) 

-0.497296223 

- 33.522088710 

5) 

-1 .296577540 

+ 32.512660280 

6) 

-1.296577540 

- 32.512660280 

7) 

-1.623612400 

•f 30.0 

8) 

-13.265108300 

+ 30.0 

9) 

-13,290158100 

+ 30.0 


3.11 COKLirSICIfS 

In this Chapter, the theoiy of aggregation due to 
Aoki is ^plied to reduce the model and compute suboptimal 
controllers for the tv^j-area 9th order 13*0 system, fwo 
oases of model reduction by aggregation are considered, 

6th order and 4th order. She suboptimal control by 6th 
order aggregation is accurate for all practical purposes 
as could be seen from figures 3.1, 3.2 and 3.3. Ihis faot 
is confirmod by a) tte increase la tte performance index 
figure which is only 1*3^ of the optimal one artd b) tte 
cloccd loop eigenvalues which are more stable than those 
obtidnod by 4th order aggregation. In tte case of 4th 
order, tte error in response is considerable; also the 
increase in tte performance index figure is 8.6% 
optimal one. 
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The suboptimal control by aggregation considerably 
reduces the computationalL effort and: tine . This could be seen 
from the fact that the matrix Eiccati eq,uation of order 
45 v/ould have to be solved if optimal control were to be 
determined, whereas in the case of 6th order aggregation, 
the order of matrix Eiccati equation to be solved is only 
21 . If tte results obtained v/ith 4th order can be consi- 
dered to be tolerable for all practical purposes, then the 
order of matrix Eiccati equation to be solved is only 10. 

In the above method correct knov/ledge of the eigen- 
vectors of the original system is essential. However, 

g 

Aoki has given a method for weakly coupled systems which 
requires only an approximate knowledge of the eigenvectors. 
This method involves an iterative process for computing 
suboptimal control. 

The above results shov^ that the accuracy of repre- 
sentation increases with the order of tin aggregc^ted model. 
Thus a tradeoff lias to be raade betv/een accuracy of repre- 
sentation and saving in computational effort and the system 
reduced to the appropriate order. 



GHAPCCEE 4 


SUBOZDUfiAL EEGUIATIOIT OP A ITOIELIMAE LOAD PEEaUEHCY 

OONLEOL SYSTEM 


In LEG systems the tieline povYer is a function of 
the sine of the angular difference of the Toltage vectors 
at either end of the tieline. Under small load disturban- 
ces .this tieline power becomes a function of the angular 
difference itself and the system can be treated as linear. 
However, under large load disturbances it has to be treated 
as nonlinear having the sine teim as the nonlinearity. Such 
LEG systems will henceforth be termed as nonlinear LEG 
systems having tieline power nonlinearity. 

In Ghapters 2 and 3 methods are presented for the 
optimal and suboptimal regulation of an LEO system. The 
analysis t lie re in is based on the assumption that load 
disturbances are small and thus linear analysis and design 
methods are ;iustified. In practical systems, load distur- 
bances are not necessarily confined to be small; and any 
analysis made should also cover large disturbances. As 
explained above , the tieline power nonlinearity is an 
important nonlinearity to be considered under such large,, 
disturbances. 
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Whereas well defined, closed form methods are 
available for the optimal and suboptimal regulation of 
linear dynamical systems, such closed form solutions have 
not been dealt with in literature for the regulation of 

7 

nonlinear dynamical systems. However, Lukes has given 
a method for the optimal regulation of autonomous nonlinear 
dynamical systems ?/hich have a linear term in control. The 
hypothesis used thei-ein is that mth the linear part of the 
system and the linear control tern, a stable closed loop 
system can be constructed by employing linear optimisation 
methods. As LKJ systems fit into the above class of 
nonlinear dynamical systems, Lukes' method is applicable 
to- the same. Section 4.2 describes Liikes' method in 
brief along v/ith statements of* some important theorems 
v/hich contain the main results, 

Lukes' method suffers from the disadvantage that 
the amount of computation increases considerably, as the 
order of tiie system goes up. Only systems of small order 
(say 6 or less) can be treated without any computational 
difficulties because apart from the usual solution of 
matrix Siccati equation, the problem reduces to the 
solution of a system of linear algebraic equations and 
the order oi the system of equations increases steeply 
as the order of the system goes up. This is explained 



jpurther in Section 4.3. Hence it is imperatiTe that some 
method of reduction he adopted for first reducing the model 
of the multiarea LHC system before Lukes' method is applied 
to the same. Section 4.4 gives a method of reduction of the 
nonlinear LHC systems, talcing advantage of the sine term in 
the tieline power nonlinearity occirrring therein. In 
Section 4,5 the method is applied to a 4th order example 
which could be considered as a single area LHC system tied 
to an infinite system . Two types of study are made on 
this 4th order example; (i) Lulces' method for determination 
of nonlinear regulation is applied to the system as such; 
(ii) the system is reduced to a 3rd order nonlinear system; 
the optimal controller determined by Lukes' method for this 
reduced system is used as a suboptimal controller for the 
original system. Thus this example serves to evaluate the 
closeness of the suboptimal response to the optimal one as 
well as the degradation in performance using the suboptimal 
controller >vLth respect to the optimal performance. 

4.2 OPTBIAL REGULATIOH OH NONLIIISM LWAillGAL SYSTEIvIS 

A. Hormulation of the Problem : 

The problem is formulated in terms of a control 
system equation in 

X = H(x, u) 

and a performance integral 


(4.1) 
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¥e seek an r-ditaeiisional yector feedback control 
function of tiie state x, u = u(x), which minimizes the 
integral in (4t2) subject to the constraint given in (4.1). 
Equation (4.1) is an autonomous differential equation where 
G(x,u) can be represented in the form 

li) = A X + B u + ^(x, u) (4.3) 

n) = S + 2x,E u + u.R u + g(x,u) (4.4) 


Here the inner product notation 
h 

x-z = yic 


and the partial differential notation f^(x,u) 
are used . 


3 f (x,u 
- - 3 U '"( 


B . De fi nit ions ; 

(i) A real matrix is called a stability matrix if all 
of its eigenvalues iiave negative real parts. 

(ii) The pair of matrices (A,B) as well as the control 
system defined by E(x,u) in (4.3) is stabilizable if there 
exists a real matrix ■ for which (A + B K) is a stability 
matrix. 

This is e^QDlained as follov® . We consider tlxj class 
of feedback controls which are of the form 

u = u(x) = K X + ^(x) (4.5) 

where h(x) denotes the higher order terms. The real matrices 
K are always selected so that u(x) stabilizes the system 
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described in (4.1)5 that is, we demand that in 

= (A + 3 K) X + 3 h(x) + f(x,u(x)) (4.6) 


(A + 3 S) should be a stability matrix . 

'/ \ w 

(iii) A C stabilizing feedback control 
i^(x) = Kjt + ^(x) is optimal for the process in (4.1) 
with respect" to the performance integral given in (4.2) if 
for every stabilizing feedback control u(x) = K x + h(x), 
thei^ exists a neighbourhood of tl:© origin in which 



(4.7) 


0, 3asic Assumptions : 

(i) It is assumed that 3(x,u) and G-(x,u) are defined 
in some neighbourhood of the origin in and can be re- 

presented in the form of (4.3) and (4.4) respectively, 

( r i ) In ( 4 . 4 ) we as sum e 




i_ -I 

(iil) Ihie fundamental hypothesis is tiiat P(x,u)is 
stabilizable , 

(iv) It is assumed that J'(x,u) and G(x,u) are real 
analytic about the origin in By this v/e understand 

the teims f(x,u) and £(x,u) in (4.3) and (4.4) to be real 
convergent power series about the origin beginning with 
second and higher order teims in (x,u) respectively. In 
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this situation we admit every h(x) in (4.5) given by real 

power series converging about the origin and beginning v/ith 

w 

second order terms. The feedback controls are called 0 
stabilizing: controls . 

D. Main Theorem and Theorem Oonccrning Truncated System 
The main results are given in the main theorem 

be low : 

Main Theorem ; 

Tor the c'^'^ stabilizable control process in 

« 

z = T'(z,u) = A z + B u + f(z,u) (4.8) 

with the performance integral 
00 

J'(Sojli) = / ^(xju)^^ 

o 

\ 

• + £(xj]i) d't 

(4.9) 

there ezists an optimal G stabilizing feedback control . 
The optimal control solves the functional equation: 

G-u^SojII^CXq)) = 0 (4.10) 

for all z^ near the origin and is unique in that 

(i) liit- is the unique 0^'^ solution to (4.10). 

(ii) Ujt- is the unique C stabilizing feedback control, 

(iii) Ujj. synthesizes tlie unique optimal open- £ d op control. 
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thermo re , l^(x) = 'K^ x + l]#(x) and 

Sq + v/herc the lov/est order terms 

are given by matrices I* and P*>- .0 depending upon only 
A, B, Qj E and 21. 


Ihe truncaxed or linear part of the djaiamic equation 
as well as the quadratic part of the performance integral 
are then formulated into a linear control problem. Ihe 

following is the statement of the theorem concerning the 
truncated system so formulated; 

Ihoorem concerninis: the Truncated System ; 

Eor the special case of the Main Iheorem in wliich 


X = A X + B u 


and 


00 




dt 


the optimal control is U)t(x) = x, v/hexe E* 
Here P*>- ^ solves the matrix equation 


(4.11) 

( 4 . 12 ) 


+ P*(A-BR“'^2]^) + (A-BE-’'l'r)P* - P* (BR'’”'b^ )P* = 

.. (4.13) - 


and is the unique positive definite solution. 


X is ^ global optimal control | finally 

J(SqJ li*) = So Sq (4.14) 

The theorem for the truncated system is first proved 

7 

with the help of intexmediate Lemmas . 
extended to tiic Main Theorem. 


The proof is then 
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E. Calculation of Po?/er Series for UjtCs.) iCxjU*) 

' -f' ' ' ' ' ' ^ ' 

Since ]^(z), J*(x)=J(x,Ux-) and E^(x) = (s) ) 

are analytic at)out the origin"* they can be expandatLn' po?/er 


series as; 

= ui''^(x) + (4.15) 

J*(x) = + ... (4.16) 

E*(x) = ^ 

The lowest order terms in the abore series have 
been computed as solutions to the truncated problem. 

£ ( 4 . 18 ) 

J^^^(x) = x.E*S (4.19) 

^(x) = X (4.20) 

where A* = (A + B K* ) (4.21) 

The computation in the case of truncated system 


reduces to solving the quadratic matrix equation 

+ E(A-BE"'’l\T^) + (A-BE“’^H®)P - P(BE"'’b'-^')P = 0 

' ( 4 . 22 ) 

for P = P.j<r-> 0 and computiiag K* and A^ by the formulas 

+ 3® P,) ■ ■ (4.23) 

and that given in (4.21), 

The follov/ing procedure is followed for computing the 
remaining terms in the power series. 

In one of the Lemmas given in Lukes theory the 
following relation is proved: 
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about the origin. 

Also from the main theorem 'm have: 
^■n (dS) ) ) + Gr (z,u^(x)) = 0 


(4.24) 


(4,25) 


After substituting for f(x,U 5 f(x)), G(z,Ujt(x)) and \(x,u*(x)), 
G^(x,u*(x)), these equations can be va-itten in the fom; 


A^x + B(u^ - Eh .x) + f(x,U3(.) . 


+ x,^ + 2x.lTujj. + Uj^.R ^ + g(x,u^) = 0 (4.26) 

u^(x) = - i R"*’ (B + J^(x,u*) + 2 A + g^, 

.. ( 4 . 27 ) 


Hence 


ihx.J^(x,u^) = - j3(u^-K^a) + £(x,ii^) .J^(x,u^) 

- 2x.Hu3(.- u^.Riit. - £(x,u^) - x.<^ 

... (4.28) 

liK-(x) = - 2 ® (S + £^) . J^(2.yU#) + 2E^x+i^(x,u,(.) 

, ( 4 . 29 ) 

Substitutii-ig the power series for u^ and J(x,^) and 
then selecting the mth order terms from the former equation 
and the 1th order terms from the latter, and using the 

r^fpv T 

equality -JB ^(x) + 2R x + 2R u^ ' ^ =0 in (4.28) the 

following equations are obtained. 
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m-1 


ife2 - ' . . . ^ ' • - 

# 

Cm-1)/2 (X) (m-l)/ N 


1=2 


- S.^“^(x,u,t) (4.30) 

* fOJP'MI ss ^i^jS**** * 


(1% 

M Hi) 


1~1 


^ ' I ■■: ! ■[ ■ ^as*! 

+ ^03? 1 = 2,3, .. . 


(4.31) 


Here fil denotes tte integer part of 1 an2 the term with 

Jm/2l ■ 


in it is to he omitted for m odd 


With the help of the ahoTre two equations and starting 
with u^^^( 2 c) « K*x and J^^^(5.) » 3c«^*5. terms in the fol- 
lowing sequence can be computed consecutively! 

I3i^^(x), J^^^(x), ui^^(x), J^^^(x),ui^^(x), ... 


thereby generating the power series for J*(x) and 3 t(x). 

$h0 computation of successively higher oriier terms reduces 
to solving successively higher oiaur systems of linear 
algebraio equations. 

F. Ixtension Of {gheorems to the Bifferentiable Case 

With some additional assumptions on f(s,,ii) , g(3c,u) and 
h(£) the proofs of the above theorems are extended to the case 
v/here i‘(£,2i) and 0(3^,^,) are twice continuously differentiable 
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ar'out the ori:;;in in and .|i(x) is at least once conti- 

nuoucl.y dilTc.rcntiable about tlie origin. 

4.3 FtEDiOTici: cr DYiii.iic/jj s'’'STd::B 

In Section 4.2 it is stated that the optimal feedback 

control&rfor a nonlinear system is of a scries form and tte 

computation of successively higher ox-der terms in the series 

uitiuatcly inuluces to solving successively higher order 

0 , 7 toms of linear algebraic o filiations. Por all practical 

/ \ / \ 

.rai’pc-.rn:; it is sufficient to stop with smd 

c>'^oau.>o fchc effect of tile higher terisa is negligible. Honce 

tl» problem further reduces to tiiat of assuming a third 

degroi form for and solving for the coefficients of 

the Ssime. Section 4.5 gives the application of the method 

to a 4th order nonlinear system and thus provides an insigiit 

into th' motiiod, Ihe order of the linear system of algebraic 

r 3 1 

enustions tiiat is to be solved for determining in 

tte case of systems of order 2 to 8 is given in fable 4.1 . 

f.^ble_ 4„.2 

Order of system vs. oru:.r of sixiultaaieous equations 


Order of the 
system 


Ord-:r of simultaneous 
equations to be solved 


2 


4 

5 

6 

7 

8 


4 

10 

20 

35 

56 

84 

120 
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Table 4.1 shows that the order of linear algebraic 
equations that is to be solved increases steeply ' as tlie 
order of the system goes up. Itirther the smallest multiarea 
IPC system, viz, a tv/o-area system is consisting ofSstate 
variables, this order increasing to much more than 8 in 
three-area and higher order multi-area systems , Hence 
considerable amount of computational difficulty can be 
expected if Lukes' method for nonlinear regulation is 
applied to tlie original model of the multi-area system. 

This suggests that if the order can be reduced by some 
means^the nonlinear regulation method provides us with a 
very good suboptinal control for the original system in 
the shape of optimal control of the reduced system. 

fortunately, in the case of LfC systems, the tieline 

power nonlinearity is suited for such model reduction, for 

6 

this purpose, the method of aggregation due to Aoki is 
suitably modified and applied (Chapter 3). This is dealt 
v/ith in Section 4.4 as applied to a two-area IfC system, 

4.4 HELUCTIOII Of IIJITIAPIIA IfO SYSTBbIS AlCD HBGHLATIOH 

The dynamics of the two-area nonlinear IfC system, 
which is linear in the control u, is given by 

x = Ax+Bu + f(x) (4.32) 

where x is an 8-vector and u is a 2-vector. 

The description of the matrices A, B and also that 
of the nonlinear part f(x) in (4.32) are given in Chapter 5. 
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2 hc ,:crforr!:ancQ index to be ninj.iiiiaed is taken as 


a;- 




o 


/ (- 1 . si + 21 u) dt 


( 4 * 33 ) 


It ia required to reduce the 8 tli order cyotGiTi dss- 
ciibed abOYc to a 5 tli order model. Recalling the theory of 
ag.,:rogation presented in Chapter 3 , the transformation tliat 
i.j made 1130 of, jior this purpose, is 


r. C 


( 4 . 35 ) 


foliomne not hod 
cciapiitc C in the 
the special type 
is nanipulatoi , 


Icncnt vector and C is a 5 x 8 mtrix. Ilio 
extends the tcclmiquo of aggregation to 
cuDo of nonlinear LFC S2f stems. For this 
of vector function f(x) occurring therein 
f(x) in the two area IiPC cystein is of the 


ions 




0 

Sin 2 w (x^ -x^ ) + k2 { Cos 2 Tr (x^ -x^^ ) - 1 } 

0 

0 

0 

k2Sin2w (x^-x^) + { 0032^^ (x^-x^ )-1 } 

0 

0 


( 4 . 35 ) 


Iho second and 


be written 


f o (x) 


''I 



sixbii nonzero eleraents of ( 4 . 35 ) can also 


Sin2ir(x^-x^) + k2^Cos2^(x-|-x^)-1} 
Sin 2 'ff (x-j-x^) 4 - k^ {Gos 2 it(xi~Xg)- 1 > 

2 ’ J 


2 ff (x^-x^) 
2w (x^-x^) 


...( 4 . 36 ) 
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let 

9-J (x) 

and 



Sin2 irCx-j -x^ ) + 1^2 ( Cos2't (x-j -x^) -1} 
_ _ 

k;jSin2'^(x.j*-x^) + k^-t Oos^f (x-|-x^)-l} 


(4.37) 

(4.38) 



These scalar functions 0-j(x) and ^ 2 ^—^ plotted 
Yersus 2Tr(x-|~x^) or ^^-j2 4^ 4he usual range, say 0 to ^/2 rad. 
A typical curve obtained for 9^ (x) with k^ and kg taken as 
0,866 and 0.5 respectively, is plotted in I*igure 4.1. It 
is seen that the variation of ®-j (x) is nearly linear in 
the range taken. The plot up to a lesser angle will 
ofcourse be more approximately linear, let the maximum 
excui-sion of ^e o radians). Por the purpose 

of computing C the value of 0.| (x) read off from the plot 
given in Pigure 4.1 at an angle midway, i.e. ^ radians 
is taken to be constant for the whole range from 0 to 
o radians, let tlje values read off at ^ radians for 
0-j (x) and 02(^^ and hg respectively. Thus the 

second and sixth nonzero elements of f(x) given in (4*36) 


can be written as 


f2(x) = h^ 2i«’(x^-x^) 
fg(x) = hg 2^(x.]-x^) 


(4.39) 


Substituting this linear expression in the place of f(x) 
in ( 4 . 32 ) and then incorporating the corresponding terms 
in the matrix A, (4.32) can be written as; 

X = H X + B u 


(4.40) 




A6|2>rad 


P»g- 41 f'lot of 0 IMj2)Vs 48|2 
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Ihe matrix H of ( 4 , 40 ) is utilized for tte puipose of 
computation of the matrix 0 using the theory of aggregation. 
Thus the five rows of G consist of the eigenvectors of the 
matrix corresponding to the five dominant eigenvalues. 


Reverting back to (4.32) and multiplying the same through 
out by 0 

Ox = CAx + GBu + C f (x) ( 4.41 ) 

Denoting 

G = C B ( 4 . 42 ) 

and noting that 

z = 0 X (4.43) 

equation (4.40) can be written as 

z = CAx + Gu + G f (x) (4.44) 

Using the generalized inverse concept, x in (4.34) can be 
written as 

X = 0^'(C z (4.45) 

Thus C A X in (4*44) becomes 

0 A X = C A 0^(0 C^) \ (4.46) 

Denoting 

1 = G A C®(C (4.47) 

and using the relation (4.45) to express Vm terms in 

^ of ( 4 . 44 ) in terms of Zj (4.44) can be written as 

z=Rz+Gu + g(z) (4.48) 

where i(z) = G ^[” 0^(0 G^)^/ (4.49) 

L_ - 4 ' 


and is a five element vector. The linear terras of z in the 
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Tector function s(z) of (4.48) are incorporated into the 
matrin P and (4-,48) is finally written as 

z=Ez + Gu + g_(z) (4.50) 

Here H is a 5x5 matrijc, G is a 5x2 matrix and £(z) is a 
5-elerrient vector fimction of z. Equation (4 #50) describes 
the dynamics of the reduced (5th order) system. Ihe 
integral performance index for the reduced system is taken 


as 





% = 

CO 

/ 


L 

% 2 + u R u) dt 

(4.51) 

0 


■ - 

- . 


where 


33^-1 




qi m —1 

C Q C-'(C C^) 

(4.52) 

and I 






; 1 

ol 



H 

i 

0 

1 

1 ! 

i 

(4.53) 


Lukes' method is now applied for the regulation of 
the nonlinear 5th order control problem described by (4.50) 
and ( 4 . 51 ). Let the optimal control vector u thus deter- 
mined for tte 5th order system be given by 

Uj^(z) = z + KgCz) + K^(z) + ... (4.54) 

A A. 

where K-j is a 2 x 5 constant matrix, £ 2 ( 5 .) is a vector func- 
tion consisting of only second degree terms in z etc. 

Using the transformation z = C x and writing (4.54) in 
terms of X only, the resulting expression, say y^Cx) 
is written as 


(x) = ^ ^ + ^2 (x ) + (x ) + 


• • • 


(4.55) 
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where 0 is a 2x8 constant laatrix and K2(x), K^(x) 

are second and third order power series in x. Thus the 
control described in (4.55) becomes the suboptiiaal feedback 
control for the oi'igiiial system. 

In Section 4.5 belov/, the above theory is applied to 
a 4th oi’dar system v/hich could be considered as a single- 
area liPC system tied to an infinite system. In Chapter 5 
computations are presented for reduction of the two-area 
LPC system of 8th order first to a 5th order model and then 
to a 3 rd order model. Suboptimal controls for the original 
system are determined in both these cases. 

4.5 EXj\i.lPIE OP A 4IH CEhSil SYSIEM 

In tliis section, optimal feedback control is first 
determined for the 4th order nonlinear system described 
beloY/ by directly applying lukes' method to the same. Next, 
it is reduced to a 3rd order nonlinear model by the method 
presented in Section 4 . 4 . Ihe optimal control constructed 
for the 3rd order model is then used as a suboptimal control 
for the 4th order system, Ihe clbseress of the suboptimal 
response to that of the optimal response is examined; also 
the degradation in performance v/ith suboptimal control, is 
evaluated . 



given by 

.^2 


0 

0 

0 


1 

t^D. 


1 


0 


2H. 

0 

1 


■1 


t 


^®gv1 


2H. 

1 

0 


4th 

order 

0 ” 


x^ 

0 


X2 

1 



^tt 


X3 

1 



^gvtj 


-^4 


0 


0 

f*Ti. . : 

“ [*1*5Ccos 2 itx^ -T ,0 )+0 .SSbSinlirx^] 

Hi 

0 

0 


0 

0 





ViJ 


(4.56) 


As stated be fore, this 4th order system could be interpre' 
ted as a single area IFC system connected to an infinite 
system. Accordingly the state variables x-^ to respec- 
tively represent/ Afdt, Af, A P and X The notation 

s 

as vvell as the parameters used here are the sauB as in 
Appendix A except for ■T2_p which is equal to 0,1 . 


Equation (4.56) is written in the vector form as 
X = A X + f(x) + B u (4.57) 

In (4.56) orhLy the differential equation pertaining to 
X2 is having the nonlinearity . The only nonzero element 
of f(x) in (4.57) can be approximated by expanding Sin 2 n 
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und Co.3 2Tr;c^ u,! to oocond order teims as 

.2 


f*T 


Ip 


211^ .. 




(?. ^ i:i ) 


2 + C.866(2itx^) 

Incorporating; tlic linear part of f(x) into the system 


matrix, (4.57) is written as 

i = I X + g(xj + B u 

Vi/ix XV 


(4.58) 


0 


t^T 


■y^^(0<856)2'" 

^ ‘ 1 ’ 

0 


1 0 

f"D-{ £*_ 


2% 


and 




0 

0 

■r^T 


0 

1 


1 


m 


®iVi 


2H. 

1 

"tl 

0 


1 


0 

0 

1 


tl 

1 


gvl 


(4.59) 


C 

0 


IP (2Trx-j)‘ 

i'l 0,5 J— 


(4.60) 


Bor the application of optimal control theory to the system 
given in (4.57) the pox’formance index is taten as 


J 


re 




Q 

ol 

X 

/ 

/ -I 


i 

■ ! 


0 

] 

0 


u 






~0 

0^ 0 ( 





dt 


(4.61) 


0 

0 

0 


1 

0 

0 


0 

0 

0 


0 

0 

0 


(4.62) 


and R is a unit scalar 
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Lukes' theory presented in Section 4.2 is nov/ 
applied to this systcni. The following features are noted 
for the system under study: 

(i) The iratrix IT of (4.9) reduces here to a 4-element 


null vector. 

(ii) g:(x,u) of (4.9) becomes zero here. 

(iii) Corresponding to f(x,u) of (4.3) here y/e have 


£(x) 


0 


f^T'i 2 

0 . 5 ( 2 ^) 2 ( 

0 



0 


(4.63) 


(iv) The matrix quadratic equation described in (4.13) 
reduces in this case to 

q + P*I + P* - ?* B B* = 0 (4.64) 


The truncated S3;-stein is first taken up for study. 

It is checked up for controllability; then matri;!: Riccati 
equation (4.64) is solved f or P^ . The same is computed as 


”0.45777^10° 

>-0.62143x10“^ 

-0.18810x10° 

-0.51001x10"'* 

-0.62143x10“'^ 

0.33741x10° 

0.26169x10° 

0. 47288x1 0"”* 

-0.18810x10° 

0.26169x10° 

0.33747x10° 

0.75492x10“"' 

-0.51001x10“"^ 

0.47288x10-’* 

0.75492x10”^ 

0.18087x10“^ 


Equation (4.18) is written hero as 
ui^ ^ (x) = Ki(. X 


(4.65) 
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and 


-1 T 
K* = -R 


( 4 . 66 ) 


where K*, a 4-eleriient row vector, is computed as 

0.63751 -0.59110 -0.94365 -0.22608 

Phe system matrix 1* = (A + B Eljj.) described by (4.21) is 
computed as 

r~o.o O.lOOOOxlo”’ 0.0 0.0 

-0.32649x10*^ -0.49980x10“'^ 0.60000x10*' 0.0 

0.0 0.0 -0.33333x10*' 0.33333x10^ 

0.79688x10*' -0.12597x10^ -0.11796x10^ -0.15326x10^ 


At tliis stage a check is made to see whether the 

system obeys the fundamental hypothesis, viz, stabilizabi- 

lity, Por this purpose the eigen values of the above 

matrix - (A + B K^) are detoimined ass 

0 


1) ^ - 0.49353x10 + jO.O 


2) - 0.13345x10 + jO.O 

3) 

4) 


3) - 0.24352x10’' + j0.58322x1o’' 


0.24352x10*' - o 0.383 22x1 o’' 


Phe real parts for all the eigenvalues are in the loft 
half plane and hence the system is stabilizable , Phe 
next step is to compute the nonlinear part of the feed- 
back controller. Equation (4.30) for m=3j here simpli- 
fies to 

f '-7 \ f r\ \ fr\\ 

(4.67) 


A»x.Jp7x) = -h2>Cx).J_l2)(x) 


X 
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III (4.67) all the teims are kno'.vn except J^^^(x) which is 
assuEBd to have tlio forn 


^( 5 ) 


(s) = 


C^X^ + CgX^Xg + c^x-jX^ + c^x^x^ + C^x^x| + CgX^X^ 
+ C^X^X^^ + CgX^ + CgX^X^ + C^qX^X^ + c^^x|x^ 

+ + '' 1 3^3 + ^= 14 ^ 3^4 5 "^ 3^4 + ° 16^4 

+ C^^X^XgX^ + G^gX^XgX^ + C^gX^X^X^ + C2 qX2X^X^ 


. . ( 4 . 68 ) 


B7 performing the products as given in (4,67) and 
equating coefficients of like powers of x on either side, 
the computation reduces to thmit of the solution of a system 
of 20th order linear algebraic equations given by 


^ (4.69) 

whejre is a 20x20, matrix and n_ is a 20x1 vector. Ihe 
description of the elements of ^ and n_ are given in 
Appendi-v n , xhc elements of the matrix v are constituted 


of the elements of the matrix which is computed above; 
hence their numerical values are not given. Finally the 


elements of the solution vector _n which are the same as 
c-j to C2 q, the coefficients in the expression for J^^^(x), 


are computed as 

c-[ = ■-0.78146x10° 
C2 = 0.19132x1®° 
C3 = -0.27530x10“ 
C4 = 0. 64460x1 


c^ = 0.38962x10”^ 
cg = 0.17074x10° 
Cy = C. 26203x1 C”^ 
Cg = -0.14959x10”^ 
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= -0.25304x10' 
= -0.21228x10" 
J^-} = -0,20398x10" 
:-l 2 = -0.47398x10" 
= —0,776^1x10 
i-iy, = -0.42168x10 


= -0.84692x10" 
C-J6 = -0.61400x1 o" 
c-jy = 0.84663x10" 

Cl8 = 0.40938x10" 

= 0.24071x10" 

Con = “8* 58479x10" 


Udfr Hx) is coaputed from (4.31) with 1 = 2. The 

ih this case reduces to 


same 




(4.70) 


( 3 ) " 

liere is a 4-element column vector cons 


is ting of 


partial differentials of Hx) v/ith respect to x-j,X 2 jX^ 
and x^ respectively, as its elements. On performing this 
partial differentiation and the matri:: products in (4.70) j 
we get 


Ujl''^(x) _ d^x^ + dgX-jX^ + djX2 + d^X2X^ + dgX^ + dgX^: 


+ drjx^ + dgX^X2 + dgX^x^ + d-j^XgX^ 


(4.71) 


where the coefiicients d-j to d-jQ are calculated as 


di = -0.10672x10 


= 0.10587x10" 


^2 = -0.32753x10" 
d^ = 0.12749x10" 


Ij = 0.11513x10’ 


do = -0.25586x10 


= 0.59248x10' 


dg = -0.1 5045x1 0^^ 


0,26355x10" 


dio= 0.36549x10* 
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As stated before, tte computations are stopped \7ith 
( 5 ) f 2 ) 

and - 2]he optimum feedback control is now 

given by 

^(x) = + U)l^^(x) (4.72) 

(1 ) ( 2 ) 

v^horo u* ' {x) and u* '^(x) arc calculated as given in 
(4,65) and (4.71) respectively, fhis expression for ^(x) 
is used in tlie system equation (4.58) and the response of 
the state variable X 2 determined. Table 4.2 gives ten 
arbitrary sets of initial conditions that are taken for 
this purpose , 

^blp. 4.2 

4tli order 33’'stein - list of initial conditions. 


SI .Eo 


1 

2 

3 

4 

5 

6 

7 

8 


Initial conditions for 


X-, X2 X3 


0.0 0.0 


0 . 01 / 21 T 

0.0 

0.349/2'^ 

0.0 

0.71 5/2^ 

0.0 

0.05/2-^ 

0.025 

0.10/2 IT 

-0.005 

0.15/2 If 

0.03 

0 , 23 / 2 ^ 

0.04 

0 

• 

0 

0.001 

0.20/2^ 

-0.008 

0.30/2 It 

-0.023 


0.0 

0.0 

0.0 

0.0 

-0.02 

0 

• 

0 

1 

0.008 

O.OO^i- 

- 0.015 

-0.001 

-0.018 

0.008 

-0.02 

0.0 

0.01 

-0.008 

0.026 

0.01 


9 

10 
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2hs responses obtained vfith tne above initial 
conditions are given as continuous curves in figures 4.2 
to 4»11» In Subsection B belov/' suboptisial responses v/ill 
b© obtained iialcixig use ox a raducca siodcl, for fcne sene 
initial conditions and v;ill bo compared T/itli tloe optimal 
responses obtained as above. 


B. Suboatimal Rorrulation by Hoduction 

In Subsection A above, tiic optimal foodback control 
cf tlio 4th order ncnlineor systora v/r?.s dotcrsiinod by direc- 
tly applying Lukoe* method to the same. In tlus subsection 
the 4tli ordci’ system is reduced to a 3rd order nonlinear 
model using the theory given in Section 4.4? Lukes’ method 
is than applied to this 3rd order system to det.rmine the 
optimal feedback controller for tho same . This is used as 
a sub optimal controller for tlic original system. 

The dynttiiics of the 4th order system is remtten 
from (4.57) as 


X = A£ + 3u + tiz) (4.73) 

v/lieio n is a 4~'vector, u is the scalar control and f(x) is 
a 4-elomont vector function. Uow reduction of the above 
system to prJ order is affected using the transformation 
z = 0 X given in (4,34). 

lor the purpose of computing C, the maximum 
excursion of the angle 2trx^ is taken as, say 0.715 radian. 








Fig. 4-5 Response of Af with initial condition X)»0’05/2TI 

X2'0-026,X3=-0-02, >^=-0:01 
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0024 
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- 0*010 


- 0*018 


tjJ 
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Figure 4.1_is made use of and the value of the nonlinear 
function 0.5(Co32Tr3:., -1 ,0) + 0,866Sin2irx^ 


the same corresponding to an aiigle which is midway between 


read off from 


the range from 0 to 0.715 radian. This value is multiplied 
Ijy -to gg-t -the linear term in x-j which is then 

incorporated in the system matrix A to foam ' matrix H. Ihe 
rows of the transf oarmation matric G now constitute of the 
eigenvectors of corresponding to its tharee dominant 
eigenvalues . 


After transformation the dynamics of the 3rd order 
system becomes 

z = i z + g u + g(a) (4.74) 

where £ is a 3-vcctor, u is the scalar control as before 
and ^(z) is a 3— element vector function of z. 


F is computed as 


F = C A 0-^(0 C ) 

(4.75) 

and 


G = C 3 

(4.76) 

A.1SO 


g(5.) = C f (x) 

(4.77) 

rp 2) 

Where the -expression G (C C ) z, 

given in (4.45) is 


substituted for x. The linear teaams in the vector 
function £(z) are then incorporated in the system 
matrix F and (4.74) vvritten as 
z = F z + G. u + g(z) 


(4.78) 
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where £(z) is now comprising of only second powers of the 
elements of z. Here the matrices C, I and & are computed 
as 


0 


-0.33002 0.21296 0.71120 -0.20279 

-0.3014S -0.42007 -0.17733 0.0 

0.71632 0.24560 0.62736 0.18160 


P = 


-0.73514x10^ 
-0.3 0623x1 o'* 

0.86271x10“ 


0.2901 4x1 o'' 
-0.78371x10° 
-0.15337x10" 



0.25349x1o'' 

0.0 


0.22700x10 



-0.88118x10""' 

0.17381x10° 

-0.10861x10*' 


Ihe performance index for the reduced system is giwen by 


J = 


Qd 

/ 


0 


u) dt 


(4.79) 


where is as given in (4.52) with Q as in C4.62) and 
is computed as 


0.81562x10 

-0.79954x10° 

-0.29835x10° 


-0.79954x10° 
0.78379x1 0'' 
0.29247x10*' 


-0.29835x10° 
0.29247x1 0^ 
0.10913x10^ 


The system described by (4.78) and (4.79) is now 
solved by Lulces’ method. As in the 41iii order case, the 
truncated or linear part of the system is taken up first. 
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At this stage the controllability of the pair is 

( 1 ) 

checked up . (z ) is determined as 

^ 2 = K z (4.80) 

A 

P is the solution of the matrix quadratic equation 

+ P P + /' P - P d R”’’ P = 0 (4.81) 

At 

Here P is computed as 

0.12585x10"^ -0.76779x10° -0.66578x10° 

-0.76779x10° 0.17583x10"* 0.51315x10° 

-0.66578x10° 0.51315x10° 0.5l866xl0°_. 


E is a 3 element row yector and is computed as 

0.1 6789x1 o'* -0.78146x10° -0.51035x10^ 

( 3 ) 

As before determined by using the equation 


2.Jp^(z) = -g^^^(z).d^^^(z) 

v/he re 

P* = (P + G E) 


(4.82) 

(4.83) 


Here P* is c onputed as 


1 1 
-0.49910x10 0.48824x10 

-0.30623x10^ -0.78371x10° 

-0.37248x10"* 0.17586x10** 


0.12056x10** 

0.17381x10° 

0.72432x10" 


The stability of the reduced system is detenain^d 
by computing the eigenralues of P# given in (4.83) ns 
1) -0.51514x10° + dO.O 
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2) -0.25936x10^ +-j0.36374xlo'* 

3) -0.25936x10^ - j0.36374x1o’^ 

All lliB aboV8 6igeii\ralues ape haviiig negatiye real parts f 
and hence' the system is stabilizable . 

hoYVj in (4 *82) all the terms are known except 
Tile following expression is assxr-ned for 

^(z) = d-jZ-j + 12^2 + d^z^ + + d^z^z^ + *^6^1 ^3 

2 2 2 

+ dyz^z^ + dgZ2Z^ + dgZgZ^ + d^QZ^Z223 (4.84) 

f *5 1 

*^z is the column Tector function consisting of the 

partial differentials of J^^bz) with respect to z^ ,Z 2 
and z^,as its elements. The vector g} is that part 

of the vector gC^) which consists of only second order 
terms. By performing the products given in (4.82) and 
equating like coefficients on eitiier side, tte solution 
reduces to that of solving a set of linear ’algebraic 
equations of lOth order viz., 

Am = jE (4.8^ 

where A is a (lOxlO) matrix and g is a (10x1) vector, 
m is a (10x1 ) column vector to be determined. Its ele- 
ments are only the coefficients d-j to d^Q in the 



expression for Here is described as 


89 


o 


o ■ o 




CM 

K^ 


CM 




(M 


K\ 




CA 

+ 

(M 

CM 

+ 


VD 

CO 






o ■ o 


CM 


O 


o 


CM 


CM 

CM 


CA 

tA 

CM 

4- 

CM 

CM 

U 


^^A 

CM 


O 


CM 

^A 

%-i ^ O 


O 


tA 


O 


o 


!A 

lA 

+ 

CM 

CM 

CM 


A 

CM 

94 

CM 


CM 

94 

CM 


O O 


KM 


A 

44 O 


O 


O 


CM 

^ O 


CM 

<94 O 


O 


O 


lA 

lA 

O 94 
A 


CM 

CM 

O ^ O 

A 


< 9-1 

A 


O 


O 


CM 
. A 
94 


CM 

94 

CM 


CM 

CM 

94 

+ 


CM 


O 


CM 

9T 

A 


CM 

CM 

9M 

CM 

4- 


9m 


CM 

* 1 — 

94 

CM 


O 


CM 

94 

CM 


O 


A 

94 

CM 


A 

A 

94 

+ 


94 

CM 


O 


A 

CM 

94 


O 


A 

r-» 

94 

A 


A 

A 

94 

CM 

+ 


94 


A 

94^ 

CM 


O 


O 


A 

94 

A 


O 


O 


A 

't — 

94 


O 


CM 


94 


O 


A 

CM 

94 O 
A 


CM 

A 

<94 

CM 


CM 

9T" 

CM 


A 

CM 

O 94 
CM 


A 

T— 

94 

CM 


CM 

A 

O 94 O 
A 


O 


O 


A 

03 


94 

O 

X 

•r4 

u 

•p 

cd 

S 

Q> 

P 

P 

O 

m 

0 

H 

Q> 

P 

0 

P 

0 

U 

05 

m 

^ 1 
•r3 

P 

P 

0 

u 

0 

Xi 


I 


If 



same is computed as 


90 


o 

• 


o 

e 


o 


CO 

CM 

K\ 

• 

r— 

o 

VD 

CD 

LA 

[>- 

« 

r~ 

o 

T— 

IAt- 
CM O 
VD r- 

O {>< 

LA 

• 

•r- 

GO 

LA 

!>• 

« 

o 

o 

X 

''M-' 

CM 

CO 

CD 

« 

o 

T— 

X 

VDt- 
A O 
O T- 
CM M 

T— • 

o 


o 


o 


o 

1 


O 


O 

1 


O 


o 


o 





T~ 

o 







tA 

V- 


r— 

iTi O 





CO 

o 







CM 

o 

!>• 

o 

CO O 













MO 

t — 

T** 

'T — 

CO 't- 

o 


o 



M 

o 


o 


o 


O 

LA 

X 

LA 

A 

X 

A X 

VD 

# 


• 


« 


• 


• 


• 


• 


# 


9 

o 


o 


o 


o 


o 


o 


O 


O 


o 













1 




1 



VO 






CO 

T— 





o* 

r- 

A O 



CO 

o 






O 





ur\ 

O 

VD O 



ir\ 

v** 





CM 

\ — 





CA 


t>- ^ 



I>* 

X 

o 


o 


tr- 

lx! 





'M- 

X 

^ X 

o 


v~* 




LA 

o 


o 


'T— 

tA 

* 


• 


' • 


• 


• 


• 


• 




9 

o 


o 


o 


o 


O 


o 


o 


o 


o 









1 






1 







VD 

V- 





VD. 

r- 

CM' 

r— 



^T— 





LTV 

o 





(A 

o 

VD 

o 



CM O 





O 

\ — 






T— 


T— 



00 T- 





CM 

M 






X 

CO 

X 

o 


CO X 

o 


o 


V- 


o 


o 


!>• 

'M- 



« 


* 


# 


• 


• 


9 


9 


* 



o 


o 


o 


o 


o 


a 

1 


0 

1 


o 


o 

00- 

r— 





VD' 

T— 




r- 

CM' 






o 





CO' 

o 



a^ 

o 


o 




OJ 

r— 





LO 

^ — 



o 

-r- 

T— 

— 





H 






M 



CA 

X 


X 





o 


o 


V- 

o 


(A 

CM 

o 


o 

• 


• 


» 




♦ 


♦ 


• 


• 


9 

o 

1 


o 


o 


o 


o 


o 

1 


O 


o 


o 




T— 



vo« 

r-* 

lA- 

r~ 





VO- 

r- 




CM 

o 




o 

CO- 

o 





A 

o 




CO 

r— 



CM 

— 

LA 

T~ 





O 

T— 




00 





w 

LA 






CM 

X 


o 



o 


MO 

MO 

o 


o 


r— 


o 

* 


* 


• 


• 


• 


• 


• 


• 


9 

o 


o 


o 


0 

1 


o 


o 


o 


o 


o 

IS^. 

Y— 





VD CM 


r— 

r— 

o 






(M 

o 





VD 

o 


o 

CO 

o 






VO 







T— 

VO 

"C — 

LA 

T— 






o 






O 

M 


kJ 

t'l 


X 

o 

• 




o 

9 

IO\ 

• 

o 

9 


o 

• 


9 

CA 

« 

« 



o 

« 


o 

1 


o 


o 


o 


O 


O 


o 


o 


o 






1 










X 





o 

o 

o 







lACM 

T- O 



t—T— 

A O 
l>“ T— 






M 







— 


o 

# 


CM X 

o 


o 


CM 

^ 1 

o 


o 


o 



IN 


LA 

• 

♦ 

o 


o 


O 


• 

o 


• 

o 


o 


0 

1 


o 


o 









CA- 

VD 

o 





^ o 
^ o 


o 


ur\ 

t'r\ 

C\i 

Id 

O 


o 


CO 

T"* 

CA 


o 


o 


TT- 

CM 

LA 

9 

V— 

H 

o 

9 

o 


• 

o 

1 

0 

o 


« 

o 


« 

O 

1 


o 


O' 


o 


o 

NArvJ 





t-CJ 



00-r- 



- 



!>• 

o 





'M- 

o 



VD 







CTi 






VD 

T”* 



T— 







v— 


o 


o 


\ — 

X 

o 


VD 

tA 

X 

o 

• 


o 


o 

9 

o 


« 

o 


% 

o 


« 

o 


o 


o 


o 


o 


o' 


0.0 0.0 0.0 0.24112 0.34763 0.97647 0.35171 - 0.61246 - 0.74496 - 0.57023 

xlo"^ x10° xio"* x10^ xlO^ xlo"* xlO^ 



9T 


B'or dsscribiiig uhe vector ^ the following equalities are 
assumed, let 


^ta ^ P11°12 P12°22 Pl3‘^32 

^ih Pl2°'l2 ^22^22 + ^23^32 (4.8?) 

^tc = Pl3°l2 ■*■ P23°22 ■'■ P-33^32 


’Where ' s 'fciie elements of the transformation 
matrix C and let 

^x ~ ^rp ^ ^ 4ir^)/(4xH^) (4.88) 

Also let the matrix H = C^CO c'^) 
elements, then 



=2 = hb ■ 

*3 “ ^to 

q = ‘^12 ha + hb 

^ ^ta ^^x ^2 ^tb 

eg = 4b^ h^l irjj + 2b^(h-[-, )^f^^ 

^ ^ta ■*■ "^l^x ^1 "^13 %c 

% ^ "^^x ^12 "^13 ^tb ^to 

^9 " 2b^^^13^^ ^tb + -^’^x ^12 ^13 ^tc 

e^O = 4b^(ti^2 ^13 ^ta ^3 ^tb ^1 ^2 ^tc^ 

.. (4.89) 
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These are oomputed as 

e-j = '-0.40294x10'' eg = 0.11663x10^ 

^2 0.23146x10° =-0.10703x10^ 

^3 = 0.301 09x1 o'* = 0.31358x10'' 

= 0.87561x1c'' e^ = 0.11070x10^ 

= -0.27312x10 '' £^q= - 0.19712x10^ 

The eleraents of the solution vector which are the s^ame as 
the coefficients of the expression for J^^^(z)are computed 
as 

d-j = 0.11032x10'' dg = -0.31170x10'' 

dg = 0.14259x10° drj = 0.27200x10'' 

d^ = -0.69438x10° dg = -0.34833x10° 

d^ = -0,28688x10° d^ = -0.67608x10° 

d^ = 0.24154x10° d^Q= 0.94776x10° 

(2) 

u* (z.) is computed from (4.31) with 1 = 2. The same in 
this case reduces to 

ui^lz) = _ Ie"'' Jp^(2) (4.90) 

On performing the matrix products in (4.90) we get the 
second order control function as 

( o ) 2 2 2 

u* (z_) = e-jZ.j + ©2^.122 + ^ ®4^1^3 ®5^3 ®6^2^3 

.. (4.91) 
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Yviiere 

ei = -0.65698x10° , = , 0.1 7269x10*' 

“ --0.34848x10° = -0.1 0832x10*' 

©3 = 0.89218x10*' eg = 0.33343x10° 

'Ihe optimal feedback control for the reduced system is 
given by 

^*(£) = ^( 5 .) ^^^( 2 .) ( 4 . 92 ) 

’ " -1 

How substituting 0 x for ^ in the expressions for 
( 1 1 ( 2 .) 

TJi* ■' (^) and u* i^) of (4.92) we get an expression in 
terms of x, say w(x) , winch becomes a suboptimal control- 
ler for the original 4th order system. The same is given 
by , . 

w(x) = X + KgCx) . (4.93) 

Here' = K C is computed as 

1^.68405x10° -0.56047x10° -0.10124x10*' -0.24779x10° 

and is a 4-element row vector. 

With the same IO numbers arbitrary initial conditions 
given in Table 4.1 of Subsection A the response of the state 
variable X 2 is determined. The curves so obtained are given 
as dashed ones in Pignres 4.2 to 4.11. 

( 2 ) 

The second order performance index '(x) is 
computed for two of the initial conditions using both 
optimal and suboptimal control, and is given in Table 4.3. 
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Tablo 4.5 

41311 order system - performance indices. 

Initial Condition vdth J^(x) y/ith sub- 

optimal control optimal control 

= O.Oi'/iTr 

0^2 = 0.001 0 . 1376fe10 “^ 

= - 0.02 

= 0.0 

= 0 . 2 / 21 ^ 

Xg = -0.008 0.4(1)185x10"'^ 

X3 = o.oi 

x^ = - 0.008 

It is scon from Table 4.3 that for the first set 
of initia.l conditions, the suboptimal index '^(x) is 
higher than the optimal one by a very small percentage 
(0.36^). This shows that the degradation in performance 
with suboptimal control is negligible and that tlie re- 
duced modol correctly represents the original system. 

But for the second set of initial conditions, the sub- 
opthmal index J^^^(x) is slightly less than the optimal 
one . This can bo explained as follows: 

As can be soon from the corresponding figure 
(Biguro 4 . 10 ) the two responses arc quite close to each 
other. To correctly compare the performance index figures 
under this circumstance, it is proper to compute J(x) 
only instead of J'' .(^) in both cases. 


0.13900x10 ^ 


0.38308x10"^ 
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4.6 COl'TGLUSICITS 

111 ijliis cliapocr a method for the suboptimal regulation 
of a. iiiulticiroa LIO sysoem is presented using model reduction 
and Lukes' method for regulation of nonlinear systcras. Ihc 
methoQ. is applied to a. 4th ordorj single area IPC system to 
construct a suboptimal controller for the same, The suboptimal 
controller so obtained is compared ’>7ith the optimal one i^hich 
is determined by directly applying Lukes' method to the 
4th order system. 

It is seen from the responses dra’TO. in figures 4.2 
to 4.11 for the ten arbitrary sets of initial conditions, 
that the suboptimal response is very close to the optimal 
one. This shows that the suboptimal controller is very 
satisfactory. The order of the matrix Eiccati equation as 
well as that of the linear system of algebraic equations 
solved with the reduced model, servo to demonstrate tbs 
drastic reduction in computation accruing out of model 
reduction. 

The above results go to show that considerable 
saving in c omputationail effort can be derived without much 
sacrifice in accuracy, by applying the techniques presented 
in this chapter to the largo signal model of a multiarea 
lie system. Such a study is attempted in Chapter 5. 



CHAPTER 5 


COI/iPUTilTIOI'T OP SUBOPTIMAL GOHTROHiBES POR 
A TWO -AREA SYSTEM 

In. this chapter suboptimal feedback controllers 

are designed for the large signal model of a two-area 

IPG system using the model reduction technique as v/ell 

7 

as lukes’ method for optimal regulation presented in 
Chapter 4. Iwo kinds of reduction of the original 8th 
order system are attempted: reduction to 5th order and 
then to 33:d order* In both cases the improvement in 
3?esponse as well as that in performance over those 
obtained by applying linear suboptimal controllers are 
evaluated. Pigure 5«1 gives the block diagram of a 
single area in the large signal model of a two-area 
IPG system. 

5 .2 TWO-AEEA NOHIIEEAR IPG SYSTEIVI 

The dynamics of the two-area IPC system with a 
step load disturbance in Area 1 is described by the 
following differential equations^ 

It = ‘h 

It ‘b = - ‘b + if - 2Bf ‘^d1 

(5.2) 
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at * ^gl = ■ 

1_ 



(5.3) 

h ‘ = 

”■ I 


- ST A X , 

^gvT SVi 

(5.4) 

ft( hf2at) 

= , 

Af2 


(5.5) 

ft ‘^2 = - ■ 

f^Dp 

A f . + 

^2 + 2 H 2 

^^g2 “21^^ ^tie2 

(5.6) 

tt ‘ ^g2 = - 

1 

®t2 

‘V + 4 ‘ 

(5.7) 

ft‘^gv2= - 

^ *f> 

m “p AXq — 

Pg^H2 ^ 

4^2 ‘^SV2 

(5.8) 


A state variable representation is made for the above set 
of simultaneous differential equations, the variables to 
be considered being /Af^dt , A f^ , & , a ,/ A fgdt , A f 2 » 

aridAZ^g* Here the notatiors are the same as those 
used in Appendix A, 

In (5*2) and (5.6) the terms and AP^^g 2 

to be expressed in terms of the state variables considered. 


Ihis is accomplished as follows: 

aP^^^^ = change in tie line power of Area 1 

= Sin(A *2 + " ^12 (5.9a) 

whe re 


■" ' * -p - and A-o Ps the nominal tieline 

12 ^rT 

angle before disturbance occurs. Equation (5.9A)can be 
written as 



AP 


tie1 


"12 


Sin ^^^(CosA A^g-I) + CosA ^2 Sin ^®-{2 


.. (5.9b) 



Using "the relation 


^5^2 = 2TT(/Af^dt -/a f^dt) (5.10) 

equation (5«9b) can be written as 

^^tie1 “ ^ 

+ Oos 2Si^2w ( / Af^ dt - / Af2^i‘t) ^ (5 .t 1 ) 

Using the relation AP^^g2 = ^12 ^^tie 1 ^12 “ “^rl'^^ra 

and i’jj-j » ^*j,2 'tbfi base pov/er in Areas 1 and 2^ - — oje^ 

^^tie2 ^12 ^12 ^ *^12 - / Af2dt)-l1 

+ Oos <5*2 Sin2‘^( /Af^dt-/ Af^dt) } (5.12) 


4uX'tfC^ 


^^tie2 ^ ^1 


Substituting (5.11) and (5.12) for aP. . . and AP in 

tiel tie2 

(5.2) and (5.6) respectively and also by redefining the 

state variables in terms of their steady state values 

f*/ff 

to eliminate the term - 2 ^ d1 in(5^)the dynamics of the 
two-area UPC system can be expressed in state variable 


form as 

X = A X + B u +,f(x) 
where A = 

ro 1 0 0 


(5.13) 


f*D^ 

f^ 

r\ 


2H| 

U 

0 

_ 1 

1 

^tT 

^tl 

-1 

n 

-1 


u 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


(5.14) 


0 

-f*Bp 

-2B^ 


-1 

®gv2®2 



loo 


B = 


0 0 0 
0 0 0 


1 




0 


0 


0 0 0 


0 0 

1 


T 

gv2 


T 


Also 


^•} (s) = 2H-J- ^12 ^ ^*51^0032^ (x-j-z^)-^ + 0,866Sin2w(x-|— z^) } 

f 2 (z) = - 2^^|2 '^2^®*^ + 0.866Sin2Tr (x-j-z^) } 

.. (5.16) 

ihe other elemenis of the vector function f(z) being zero. 

The expression in the brackets[Jin both cases is the same 

and can be denoted as, say,^(z). The function , 

[_2»(zj-x^) 

is already plotted versus 2if(x^'-Z^) in Bigure 4.1. 

This fact will be made use of in Section 5.3 below. 


Bor posing the above system as a control problem, 
the performance to be minimized is taken as 

CD m m 

J = / (s 'O' S + A, S u) dt (5.17) 

o ^ - 

5.3 SUB0BTIM4L ESOUIATIOU USHG MODEI EEBUCTICIJ 

The theory for reduction of a two— area 2jKI system 
is given in Section 4.4. The next step viz., the regulation 
of the reduced (5th order) nonlinear system will be consi- 
dered in this section. The computational results in 
respect of reduction to a 5th order model are dealt with 
in Section 5,4 and in respect of reduction to 3rd order, 
in Section 5.5. 
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The starting point for the regulation of the reduced 
(5th order) system are equations (4.50) and (4.51) wiiLch are 
reproduced below for convenience: 

2=:EZ + Gu + £(z) 


TO j, 

% = / (2. % 2 + U R u)dt 

The initial condition is given by 
z(o) = 0 x(o) 

Equation (5.19) can be written as 


(5.18) 

(5.19) 

(5.20) 


00 


J. 


M 


/ 


T * T ‘ 

H. I 


e. 


2 

P E 


Ji- 


dt 


(5.21) 


Equations (5.18), (5.20) and (5.21) constitute the 
control problem for the 5th order model. It is assumed that 
the system described herein obeys the four basic assumptions 
made for the application of Inikes' iiKjthod of regulation 
(Section 4.2). hov/ v/e shall apply the theory for optical 
regulation to this control problem. 


The first step is to solve the tiuraated or linear 
part of the system and corresponding quadratic performance 
index lor obtaining a linear feedback control. Thus in 
(5.18) the linear part is written as 

i = E z + G u (5.22) 

!Tov 7 (5.22) together with (5.20) and (5.21) consti- 
tutes the control problem for the truncated system. Ths 



controllability of tbe (F, G) pair is checked up. For 
determining the optimal feedback control parameters for 
this prob lem jthe matrix Riccati equation described in 
(4.22) is solved which in this case reduces to 

F + F P* ~ P# G r”"* G^ P^(. + = 0 (5.23) 

The linear part of the optimal conticl for the nonlinear 
system thus becomes 

ij'‘^(a) = I z = -R'""' G® P^ z (5.24) 

and 

= a(o)’^ P* z(o) (5.25) 

At this stage the eigenvalues of the closed loop matrix 
F* = (F + G K) are computed and the stabilizability aspect 
of the system, checked up. The next step is to find the 
second and the higher degree terms for the feedback control 
and performance index given by 

Hs^(z) = ui^^(z) + ui^^(z) + ... (5.26) 

J*(^ = J^(z) + J^^^(z) + ... (5.27) 

•f 

As fairly accurate results can be obtained by approximating 
the controller in (5.26) as the sum of the first two terms, 
the computations can be stopped with vz) aM 

The equation corresponding to (4.30) reduces in 


this case to 
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( 2 ) 

In ( 5 « 28 ) £ , (s) = £(^) because care was taken in 

Section 4.4 to incorporate any linear terms in the P matrix, 
and thus only second degree terns are left in £(2 ) . Also 




2 P 112 -| + 2 p^ 2^2 + 2 p^ 3 Z^ + + 2 p^ 5 Z 3 

2P22^2 + 2p^2^1 + 2P23Z3 + 2P24Z4 + 2P25Z5 

2P33Z3 + 2p^3Z^ + 2P23Z2 + 2p3^z^ + 2P33Z3 

2P4424 + 2p^^z^ + 2P24Z2 + 2 ^ 34^3 + ^P 45^5 

2P33Z3 + 2p^3Z^ + 2P25Z2 + 2P33Z3 + 2p^3Z^ 


(5.29) 


where p^^ ‘ s are the elements of the solution matrix P^e- 
of (5.23) 'and 


= (P + 0 K) (5.30) 

Ihe only unknown term is which is a yector valued 

function having the partial derivatives of J4^^(z) with 
respect to to Z3 respectively as its elements. Hence 

(3 ) 

^ is assumed to take the follo?/ing third degree fom: 

t(3)/\ 3 2 2 2 2 2 

v^j — d^z^ + d2Zi-j22 + + ^4^1 ^4 "b ^2 

+ d^^z| + dgZ^z^ + d^2^z| + d^Qzl + *^11^2^3 

+ ii2^2^4 ^13^2^5 *^ 14 ^ 2^3 ‘^15^2^4 ■** *^ 16 ^ 2^5 

+ + di 8 ^ 3^4 *^ 0 ^ 3^4 ^ 21 ^ 3^5 

+ ^22^ + ‘^23^4^5 ■*■ *^24^4^5 ^25^5 ^6^1^2"^3 

+ d2y2-jZ2Z^ + d28Z'jZ2^5 ^29^1 ^3^4 %0^1^3^5 

+ d3-|Z^z^z^ + + d 33 Z 2 Z 3 Z 5 + ^34222425 

+ %5^3V5 
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Substituting lor from (5.29), E* from 

(5 *30) and , (^) oouained from (5.31) in1:o (5.^8) and 
equating xilce povvqns of _z on both sides, the equation ulti- 
raately reduces to the solution of simultaneous equations in 
35 unlmown?. Ihe same is written as 

A a = G (5.32) 

Ii6PG ^ is 3. (35xj)5) c oGf fic iguIj inBipixj ^ is b, 35 6101110111 # 
vector and ^ is the right hand vector 7 /hich is known. The 

A 

description for the matrix A an(3 that for the vector e 
are given in Appendix E . After determining the solution 
vector m of (5>32) the second degree control function u^^^(^) 
is computed* Here u'' "^(z_) is a vector valued function 
consisting of two elements. With 1=2, (4.31) in this case 
reduces to 

ui^hz) = - I G® Jp^z) (5.33) 

All the quantities on the right hand side of (5.33) 
are known, llie optimal feedback control for the 5th order 
system is now given as -the sum of the right hand expressions 
in (5.24) and (5.33) 

u*(z) = u^^z) + 4^^^) + ••• (5.34) 

or " 

Ujt(z) = -E*"^ C-^ P z - Ie"'* G^ Jp^(z) (5.35) 

Substituting 0 x for z in the right hand side of (5.34) 
resulting expression, consisting of first and second degree 
terms of the elememts of the vector x, becomes a subopt ii»l 
control for the original 8th order system. 



5.4 EEDUGTIOrr TO 5TH OETiiiE - GOEIPUTiTCIONiU:! EESUITS 


The compul; elision sijar1;s with sssuitiptioii of the 
range of the excursion of tte angle = 2ir(/Af^dt-/ Af^dt) . 

This is achieved by assuming an initial condition for the 
variables / Af^dt and /Af^dt, By applying the above descri- 
bed feedback control, the angle. A 6^^ starts from the given 
nonzero initial value and ultimately reaches a zero steady 
state; and hence the range is taken as from zero to this 
initial value . Eor instance one set of initial conditions 
taken for the state variables is 


x^ (o) 


0.30 




0.01 

“ 2ir 


with the rest of the state variables of the 8th order 
system assumed to be zero. Thus the initial angular 
perturbation A 6 ^2 = 2 T(/Af^dt - /a f^dt) ;= 0.29 and the 
range of A <S-j 2 is from zero to 0.29 radian. The next step 
is the computation of the H matrix of ( 4 , 40 ), To accom- 


plish tills the values of the scalar functions 0-j(x) and 
©^(x) corresponding to 0,145 radian ( = 0 , 29 / 2 ) are 
obtained by making use of Eigure 4»1 . Using these values 


for 

(x) and 

®2(£) 

the H matrix (8x8) is 

computed 

as 

0.0 

1 .0 

0.0 

0,0 

0.0 

0.0 

0,0 

0.0 

- 3.123 

-0.050 

6.000 

0.0 

3.123 

0.0 

0 

• 

0 

0 

0 

0.0 

0.0 

- 3.333 

3.333 

0.0 

0.0 

0.0 

0.0 

0 

• 

0 

-5 .208 

0.0 • 

-12.500 

0.0 

0.0 

0.0 

0 

• 

0 

0.0 

0 

• 

0 

0.0 

0,0 

0.0 

1 .0 

0 

• 

0 

0.0 

3.123 

0.0 

0.0 

0.0 

- 3.123 

- 0.05 

6.000 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-3.333 

3.333 

: 0.0 

0.0 

0 

• 

0 

0.0 

0,0 

-5 .208 

0.0 - 

12.500 
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The eigenvalues of the matrix H are computed as 


1) 

'1 

-0.15883x10 

+ 

O'O.O 

2) 

-0.24247x10''^^ 


DO.O 

3) 

-0.13290x1 0^ 

+ 

do.o 

4) 

-0.13266x10^ 


DO.O 

5) 

-0.12966x10"* . 

+ 

30.25127x10"* 

6) 

-0.51441x10° 

+ 

30.34764x10"* 

7) 

-0.12966x10"* 

- 

30.25127x10^ 

8) 

-0.51441x10° 

- 

3 0.34764x1 o'* 


The five dominant eigenvalues bearing numbers 2), 6), 

8), 5), 7) in the above eigenvalues are considered and 
the corresponding eigenvectors are computed, finally the 
5oc8 matrix C the rovfs of which are constituted of these 
eigenvectors, is computed as 


0.544 

0.213 

0.384 

0.102 

0.544 

0.213 

0.384 

p.107 

0.435 

-0,111 

-0.369 

-0.103 

-0.435 

0.111 

0.369 

0.103 

0.265 

0.264 

0,107 

0.0 

-0.265 

-0.264 

-0,107 

0.0 

0.0 

0.144 

0,588 

0.175 

0.0 

0.144 

0.588 

0.175 

0.0 

-0.291 

-0.132 

0.0 

0.0 

-0.291 

-0.152 

OvO 


The matrices 1 (5x5) and G (5x2) of (4.50) and the matrix 
% (5x5) of (4.52) are computed as 


F=. 


0.0 0.0 0.0 0.0 

0.0 -0.48089x10° 0.55215x10"* 0,0 


0.0 

0.0 


0,0 


-0.35559x10 


1 


0.0 


0.0 


-0.62135x10° 0.0 


0.0 


-0.12966x10^ 

-0.25127x10^ 


0.0 


0.0 


0.0 

0.25127x10^ 

-0.12966x10^ 


0.0 
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G 


0. 12802x10"^ -0.12832x10"^ 0,0 0,21884x10^ 0, 

p .1 2802x10"^ 0.1 2832x1 o'* 0.0 0.21884x10"^ 0. 


0 




% 


0,393x10 

0.0 

0.0 


-0.1 07x1 o' 
-0. 524x10* 


-2 


0.0 

0.0 

-0.107x10"^ 

-0.525x10' 

0.1 65x10° 

-0.552x10° 

0.0 

0.0 

■0.552x1 0° 

0.185x10^ 

0.0 

0.0 

0.0 

0.0 

0.293x10° 

0.143x10 

0.0 

0.0 

0.143x10^ 

0.700x10 


The next stop is determination of the feedback controller 
consisting of first degree and second degree terms in the 
elements of the vector z. 

The truncated control problem described by (5.22) aru3 
(5.21) is then solved. At this stage the controllability of 
the (T',G) pair is checked up. Uext, taking R as 2x2 identity 

A 

matrix, the solution matrix P* of (5.23) is computed as 


5x5 matrix 




” 0.2 29x1 O*"^ 

0.0 0.0 

0.507x10“^ 

-0.151x10**' 

o 

• 

o 

0.509x10° -0.165x10° 

0.0 

0.0 

0.0 

-0.165x10° 0,478x10° 

0.0 

0.0 

0.5 07x1 0**^ 

0.0 0.0 

0.335x10° 

-0.328x10° 

-0,151x10 

0.0 0.0 

-0,328x10° 

0.1 66x10^ 

The linear feedback matrix E (2x5) 

of (5.24) is 

computed £ 

-0.140x10”^ 

0.653x10° -0.212x10° 

-0.733x10° 

0.720x10° 

-0.140x10"^ 

-0.653x10° 0.212x10° 

-0.733x10° 

0.720x10° 
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The matrix 

of (5.30) 

is computed 

as 


-0.359x10 

0.0 

0.0 

-0.1 88x1 o'* 

0.1 84x1 o'* 

0.0 

-0,21. 6x1 o’* 

0.407x1 o’* 

0.0 

0.0 

0.0 . 

-0.356x10'' 

-0.621x10° 

0.0 

0.0 

-0,614x1 0 ^ 

0.0 

0.0 

- 0 . 450 x 10 ’' 

0,567x1 o’* 

0.0 

0.0 

0.0 

- 0 . 251 x 10 ’' 

-0.130x10^ 


Tho st abilizability aspect of the system is nov? 
checke.d up by computing the eigenvalues of P* as 

1 ) -0.14294x10"^ + jO.O 

2) -0.2901 7x1 o’* + ■30.34142x10’' 

■ 3) - 0 . 29017 x 10 ’' - 30 . 34142 x 10 ’' 

^ i "t 

4) -0.13885x16 ' + 30.37243x10 

5) -0,1 3885x1 o’* - 30.37243x10^ 

As all the eigenvalues are having negative real parts, the 
5th order system is stabilizable . 

The coefficient matrix A and tiie right harod vector £ 
of ( 5 * 32 ) are computed and given in Appendix E. The ele- 
ments of the solution vector m (35x1 ) are computed as 

^10 ~ ”0.44360x10° 
m-ii = 0.56631x10"’' 

m-j 2 = -0,12688x1 o’* 
in-jY = 0.9.6424x10° 

foe rest of the elements being zero. 
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as well as Ip^(z) are now completely detormined. 

The 2-element vector feedback function is finally 

determined as 

O.53455Z2 + O.OlZeiz^z^ - 1 . 6281 z 2 Z^ 

2.2422z2^ - O.07267Z2Z3 + 1.6281z2Z^ 

Now the optimal control vector function %(z) is given 
as 

li*(23 = 
where 

^ (z ) = K z 

as given in (5»24) and is computed as above. By 

substituting £ = C x in the above optimal control vector, 
the same vvill become the suboptimal control vector say 
y*,(x) (2x1) for the 8th order IBC system. 

The responses of the state variables Lt-^ and A f 2 
of the 8th order LBC system with y^Cx) as a feedback 
control function are determined and are given as continuous 
curves in Bigure 5.3 • i’or comparison, the system is lineari- 
zed about. the operating point, the linear suboptimal control 
function determined by 5th order aggregation in the same 
lines as detailed in Chapter 3 and applied to this nonlinear 
system. The same is given as dashed curve in Bigure 5.3. 

Two other initial conditions are taken, (a) x^= 0,06/2*, 

X3 = 0,00l/2Tr, rest of the states zero and (b) x^=1 ,099/2*, 

X3 = 0,O174/2ir, rest of tte states zero? and. the suboptimal 
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responses of the state variables a f-| and & are deter- 
mined as shown in Figures 5.2 and 5.4 respectively. Also 
the responses obtained by applying the suboptimal control 
assuming tloe system to be linear are shown as dashed curves 
in Figures 5.2 and 5.4 respectively. Ihe perfomance index 
figures J (x) obtained with the suboptimal nonlinear 
control and suboptimal linear control and determined for 
the three initial conditions considered above are given 
in Table 5.1. 

Table 5.1 


Seduction to 5th otder - performance indices. 


Initial condition 

(2 ) 

' (x) with nonlinear 
suboptimal control 

J ^ ^ (x) with 
linear sub- 
optimal control 

x-^ = 0.06/2 IT 

= 0.001/2’r 

0.59226x10"^ 

0.95747x10"^ 

rest zero 



x-j = 0.30/2 
x^ = 0.01/2 If 

0.13136x10"^ 

0,2l120:d0"^ 

rest zero 



= 1 ,099/2'>r 
x^ = 0.0174/2 ir 

0.10183x10""' 

0.17237x10*"^ 

rest zero 




5,5 REDUCTION TO 3BD OEFEE - CQMPIJT.'2I?ICMjaj HESSIffS 

The computational results for the reduction of the 
two-area IFC system to a 32'd order model and for the 
corresponding suboptimal nonlinear regulation^ are given 
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in this section. The matrix xH {see (4.40)) remains the 


same as computed in Section 5.4. The three dominant 
eigenvalues out oi the eight ones computed are bearing 
numbers 2), 6) and. 8) therein. The rows of the matrix 
0 (3x8) in this case consist of the eigenvectors corres- 
ponding to these eigenvalues. Thus the three rov/s of the 
matrix G are computed as the first three rows of the 
C (5x8) matrix computed in Section 5.4. 

The matrices P (3x3) and g(3x2) corresponding to 
(4.50) and the matrix (3x3) corresponding to ( 4 . 52 ) arc- 
comp utod as 


0.0 


P = 


0.0 


0.0 


0.0 0.0 

-0.48089x10° 0.35215x10*^ 

■4 

-0.35559x10 -0.62135x10° 


1 1 1 

0.12802x10 -0.12832x10 0.0 | 

1 1 I 

0.12802x10 0.12832x10 O.Oj 



-1 

0.91212x10- 0.0 

0.0 0.16488x10° 

0.0 -0.55199x10° 


0.0 

-0.55199x10° 
0.18480x10*' ■ 


The truncated 3i‘d order control problem corresponding to 
(5.22) and (5.21) is solved. The controllabilit;!/ of tr© 
pair (P, G) is checked up. The matrix E is a 2x2 identity 
matrix here also. The solution matrix P* of the 33:*d order 
equation corresponding to (5.23) is computed as the 



matrix 



0.16681 

0.0 

0.0 

0.0 

0.50855 

-0.16534 

0 .0 

-0.16534 

0.47753 


Hiie linear feedback matrix E (2x3) of (5.24-) and the 


matrix E* of 


E= 


- 0.21356 

-0.21356 


(5.30) are 

0.65256 

- 0.65256 


computed as 

-0.21216 

0.21216 


E* = 


■0.54681x10 

0.0 

0.0 

0.0 

-0.21556x10”’ 

0.40660x10 

0.0 

-0.35559x10”’ 

-0.62135x10' 


Tile stabilizability of the system is checked up by comput- 
ing the eigenvalues of E* as 

1) -0.5468fx10° + jO.O 

2) -0.55112x10° + j0.256l7xlo"' 

3) -0.55112x10° - 3 0.25617x10”’ 


As all the eigenvalues are having negative real parts, tlxj 
system is stabilizable . 

The linear part of the controller is given as 
u*'’\z) = E £ and the closed loop system incorporating_ 
this control is given as z = P* z. 

The equation corresponding to(4.30) reduces in this 


case also to: 
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Here H* is coLipuied as above. = g(z) as 

beiore and cofisiGis of only second degree terms in the 
elements of z , 


( 2 ) 


Ja .(2) = 


2 Pl1=^1 + 2P12^2 + 

2P22^2 

2P3322 + + 2P23Z2 


(5.37) 


where p^^^ ' s are the elements of the matrix P* computed 
as above, Ihe only unknown term is and corres- 

f'Z\ 

ponding to this, is assumed as 

J^ 5 \z) = lizf + 12^2 + 13^3 + 14^1^2 + Vi ^2 Vi ^3 

+ 17^1^3 + 18^2^3 ■*■ ^9^2^3 ^10^1 V3 

.. (5.38) 

- 2 2 2 - 

31-jZ-j + 2l4Z^Z2 + 1^Z2 + 21gZ^Z3 + 1^Z3 + l|o^2^3 

312^2 ^ 4^1 215^1^2 ^^^ 2^3 V3 ^ 10 ^ 1^3 

3I3Z3 + IgZ^ + 2'lrjZ^z-^ + ^^2 '*’ ^^9^2^3 ^10^1^^ 

.. (5.59) 

Substituting the above terms in (5.36) and equating 
like powers of the problem reduces to the solution of 
simultaneous equations in 10 unknoinas as given by 

A m = c (5 .40) 


Jp^(z) = 
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Here the description of the 10x10 matrix h is the same 
cis txiat in (4 •86). HoF/ever, f . are the elements of 
me-itrix ooniputod in this section as above. The elements 
^1’ '■*’ '^■JC ^ight hand vector £ are described 

below. Let 

(0.5f* %2 

(0.5f* 

q_2 = (0^2 - Cgg) 

(0.5f* ^12 

= __ _ - G^g) 

vfhere c-.*s arc the elements of the C matrix. Also let 
J 

g-| = t-j -j , S2 = ^"[2 ** *^52 ^ ” "^15 ** *^53 

y/liere arc the elGnients of the (8x3) matrix =0^(00^*) 

Then 

2 

— 2s^(q.|P-j>j + <12^12 

2 

^2 ^ 2s2(q-jP-12 + a2P22 ^ ^^23^ 

2 

= 2s2(q^p^^ + 0.2^23 *^3^33^ 

How donoting 

y-j = q-jP-j-j + 12^12 ■*■ ^^13 

^2 ~ ^1^12 *^2^22 %^23 
and 

y^ = + 12^23 ‘^3^33 2 -kxi^ 

2 

^4 = 4^1®2^1 ^2 



2 

eg = + 2s^y^ 

2 

e^ = 2s^y-| + 4s^s^y^ 

Eg = 4s2S3y2 + 2s2y3 
Eg = 2s|y2 + 4s2S3y3 

^10 ~ 4-(®2®3^1 s-]S3y2 + s^Sgy^) 


. .. (5.42) 

Here ijiiG cocff'icieiLt matrix A is computed as given in 
(5.41). The elomonts of tlie right hand vector £ of (5.40) 
are computed , as 

C2 = 0. 26673x10*' 4 £3 = -0.69565x10*' , 

^8 = 0.33333x10*' , Eg = -0.55140x10*' 


the rest of the elements being zero. Also the elements 
of the solution vector, m, are computed as 

EI2 = "0.52931x10° , m3 = 0.10679x10*' , 

mg = 0.21249x10° , mg = -0.12213x10*' 


the rest of the elements being zero. Thus the expression 

for ^(£) given in (5.39) is now determined fully. 

(2 '' 

T’inally the expression for ujt ^(z) c or i-t; spending to (4.31) 
with 1=2, is written as 


i4^^(z) = - Jp^(z) 


(5.43) 
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and is determined as 

r-1 .018822 - 0.783622 + Q.2^21z^z~ 

2 P 

1 . 018822 '+ 0 . 78362 ^ - 0.2721z^z^ 

As before, the optimal control vector 1 ^( 2 ) (2x1) for the 
reduced system is given as 

li^(2) = + iij[2)(2) (5.45) 

(1)^ (2)* 

where a and u^^(s) are computed as above. By 

substituting z = 0 x in the expression for the abo-ve 
vector function, the resulting vector function say y (x) 
(2x1) becomes a suboptimal feedback control for the 
original 8th older system. 

The suboptimal response of the two state variables 
Af-j and Af 2 of the 8th order system with each of the throe 

initial conditions considered in Section 5.4 are given as 

continuous curves in Figures 5.5 to 5.7. For the same 
initial conditions, the responses obtained with linear 
suboptimal control determined by linearising the 8th order 
system about the operating point and using the 3rd order 
aggregation arc given as dashed curves in Figures 5.5 to 
5.7. 

As in the 5th order case, the performance index 
figures J^^^(x) v/ith the above two kinds of suboptimal 
control are computed and given in Table 5*2. 



Suboptimal responses of two-area system 

ucfion 8 to 3) With initial condition 
I00l/21[,rest zero ,"on,x,-u-ub/ 2 ii 


Wit h nonlinear control 
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Table 5»2 

Reduotion to 3rd order - performance indices. 


Initial condition 

(2 ) 

■^(x) with suboptimal 
nonlinear control 

J^^^(x) with 
suboptimal 
linear control 

oi 

li 1( 

0 . 06/2v 
0.001/2 If 

0.74985x10“"^ 

0.66922x10“"’' 

rest 

zero 



X-j = 
X5 = 

0 , 30 / 2 ^ 

0 . 0 l/ 2 ir 

0.16964x10“^ 

0.15166x10“^ 

rest 

zero 




1.099/2 




0 . 0174/2 ^ 

0.15415x10“”’ 

0.13659x10“”’ 

rest 

zero 




5.6 COICLUSIOITS 

A subopt imal controller is determined for the large 
signal model of the two-area system by reducing the same to 
a 5th order model and tiien applying luleos' nethod for opti- 
mal regulation to the reduced model. Secondly the system 
is reduced to a 3rd order model and suboptimal controllei 
constructed as above. Per comparison line am suboptimal 
controllers are also determined by linearizirig the system 

about the operating point and reducing this linear model 

' 6 

by maJcLng use of the method of aggregation due to x^oki . 
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Prom Piguros 5.2 to 5.4 it is seen that the sub- 
optimal responses obtained by the method presented are 
superior to the linc^'r sub optimal controllers, Highcr 
Yalues for performance incidcs obt>-*incd with linear con- 
trollers (fable 5.1) also confirm this fact. 

Prom Pigures 5.5 to 5.7 and also fable 5.2 it is 
seen that no adTantage is gained by model reduction to 
3rd order and applying Lukes’ method to tl® same, fhis 
not only confirms the fact tliat accuracy of representa- 
tion increases with the order of the reduced model but 
also shows that for the case of two-area LPG system 
considered horo , the appropriate model for reduction 
is that of order 5. 


Optimal coartrol of nonlinear dynamicfJ. systeras had 
been hitherto accomplished by the well knovm method of 
linearization of tlie system at every step* fhis involves 
considerablo amount of computationj increover in many 
cases, the results obt:.ined may not be acc’urafe. Luo to 
the reasons stated above, the method for subop timal 
control presented in this chapter is quite promising for 
systems having trigonoraetrical sine type of nonlinearity . 



CHilPTER 6 


DWJMQr EPEECSS OB' EXCITATI02I COETROE IN A 
LOAD IREQUENCY COI:TTROL SYSEEM 


In tlie previous chapters the 'IPG system is studied 
with the assumption that there is no interaction bet'ween 
the Megawatt frequency and Ilegavar voltage control loops 
or equivalently assuming that the excitation system has 
no effect on the dynamics of the IPG system. Ihis assum- 
ption is permissible only when ti^e speed of the excita- 
tion system is much faster than the IPG loopj but in 
practical systems this is not true. Ihe importance of 
damping due to excitation control on the stability of 
synclironous maciiines lias been recently realized ^ and 
the same has been studied both by simulation and by conduc- 
ting actual field tests. lurick^^, is probably the first 
to investigate the damping effects of voltage control on 
IPG systems. He studied the effect of area voltage on 
the tieline power flow for a two-area IPG system; the 
assumptions made in his study are; (i) the area voltage 
is available as a control variable which can be manipula- 
ted for obtaining optimal response; (ii) that the area 
voltage perturbation does not have ai^^ effect on the area 


load 
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Assumption (i) made by Durick is not tinxe because 
voltages are controlled in actual systems by means of voltage 
regulator and exciter loops consisting of the amplifier, 
exciter, stabiliser am the field time constant. Assumption 
(ii) again is not coirrect because it amounts to neglecting 
an important interaction in 15’C systems. 

In this chapter a two-area EPG system is studied 
without making the above assumptions. 

Section 6,2 gives a brief description of the system 
studied by Duiick. In Sectioias 6.3 to 6.5j a more realistic 
model is treated by including the exciter loop as detailed 
above; also Section 6,6 gives the computational results 
for this model, 

6.2 OPTBIillY VOIIAGS DMEBD IT/O-AEEA SYSTEM** 

The block diagram of the two-area HO system ?fitli 
voltage perturbations A | V-jj and A |V 2 | as additional inputs 
is given in Eigure 6,1 . The system differential equations 
for this system remain the sanB as given in Appendix A 
except those for Af-j and A f 2 ; moreover tie st&te variable 
/AP^^^^dt has not been considered in the present model. 
Considering a load disturbance A in Area 1 only, the 
differential equations in respect of * f-j and * ±2 ^ 

follows . 

d ( * f "I ) f f * A f* f* 

dt : 2H| - 1 ^ 2Sf. g1 2E| -^dl 2Sf . tie1 

( 6 , 1 ) 


« « 




Figi 6-1 Block diagrann of a two-area system with voltsge 
oerturbations A|V,| and AlVjl added cs extra inputs. 
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dCAf^) Dg 

dt ■ aHg ^ ^2 + 2 H 2 ^ ^g 2 *■ 2 H 2 ^^tie 2 

The e::cpression for assumed positive in the direc- 

tion 1 to 2 is 


Yi IV 2 


Sin( 6 ^ - 62 ) 


-tlel 

If ve assume that the four quantities 6 -j> S 2 ^ 

I Y 2 j change by the amounts {"^li ^ j ''^2 

respectively, the foimula forAp_^^^^ becomes 


(6.3) 

1 t and 


‘hlel 


1 , , 31 * 4 - • 1 , 3 1 * 4 .- 4 

txelA 4. ■ ■ Jagl AY 4. _li§i_A (6^-6 ) 

M ^ 3 |Y2^ ' ^ ^ ^ 


.. ( 6 . 4 ) 

Upon perfoiming the partial differentials called for, 
using the xel&tions A5j=2ir J ^ C^*^) becoraes 


&P 


tiel 


A Y 


. . (6.5) 


where T *2 is the synchronizirig coefficient as defined in 
Appendix A, The coefficients T| and Tg are defined as: 


nr^ 

■^1 


m* 

^2 


V* 


* 


^1^12 


Yf 


^r1^T2 


3in( A -1 


Sin( 6 - 


^ 2 ) p . n . ill/p . u , V ol t 
* 

Sg ) P .U . ll'.f/p .u .volt 


( 6 . 6 ) 


(6.7) 


Substituting right hand side of (6.5) for 
APf^ei ( 6 . 1 ) tile latter becomes 
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d( A ) f* 2' 


* 

12 


dt 


- ' 2H. 


f* fyi* 


/a f^dt - 






f* 


2 a,-" * 2H7''^g1 


i.| 


Y1 


1 

f* H* 

1 

A -D 

2H| ^ ^d1 


f* tJ 

2H| 


Yr 


( 6 . 8 ) 


Hov/ using the relation 


*^tie2 "■ ^12 ^-tiel ’ 

and that in (6*5) a similar differential equation can be 
written in the case of (6.2) also. 

In (A. 5) of Appendix A the wariablesA and A 2 
are the usual LPC system control variables. lo study the 
effect ofAjv^l andAj? 2 j here the control variables 

A P^^ andAP ^2 ai'® considered to be absent. ¥ith the above 
changes tlie system dynamic equations can be m'itxen in the 
form 


x=:Ax+Bu+r aP^ 


( 6 . 10 ) 


y/here x is an eight element state vector consisting of 
the variables /a f,| dt, Af^, AP^^, A , /a fgdt, A f^, 

its elements. Also u is a 2-eleraent 


«'g2 

control vector consisting of A ly. 


ina 


B = 


0 - 0 0 0 - 0 0 


72|* Here 
T 


2H| 


2 H 2 


f*P* 

0 - 0 0 0 0 0 


2Hf 


2E2 


( 6 . 11 ) 
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A= 


0 1 0 0 


2H- 


1 


0 


0 


0 


2H-j- - - 2H-J 
0 


Vl®T 


0 




2H2 

0 

0 


0 

0 


=1 

14. 


0 


1 


0 0 


tt 


^gvl 


0 

0 

0 

0 


2H. 

0 

0 

0 


1 


0 

c 

0 

0 

0 


0 

0 

0 




2Hg, 


0 

0 


2 ^- 


2 


«2 


0 

-1 


=1 

T. 


0 

0 

0 

0 

0 


1:2 

0 


1 

^'-b2 

-1 


'gv2^2 


■gv 2 


( 6 . 12 ) 


To apply the optimal control thcorjr to tho aboTo 

system the third temrA^ in the right hand side of 

(6.10) is to be eliminated. To accomplish theis , the sans 

3 

procedure as was described in Slgerd and Tosha is used. 

That is, the states and controls hawe been redefined in teirnis 
of their steady state values. Thus tho dynamic equation 
becomes 

z = A X + B u (6.13) 

To pose this as a control problem, tlic cost function is 

taken as 

1 -7 7 

J = ^/ (x ^x + u H u) dt (6.14) 

The matrix Q is chosen such that it gives woightage for the 
frequency deviations * f-j and ^^2* same is given as 



129 



aiid. E is talcen as 

~0.0l 0 

E = (6.16) 

0 o.olj 

From Pontryagin’s minimum principle (Ap.jendix 3) 
the optimal inputs A ; j and A Y 2 I arc computed as follows: 

A Y^j = 3.575x10“^( / Af^dt -/a f^dt) + 4.699:Af^ 

+ 2.22 AP^^ + 0.279 - 4.71 A f g - 2.24APg2 

- 0.284 AXg^ ( 6 . 17 ) 

A 72 = 3.376x10“^(/Af-jdt-/Af2dt) + 4.699 Af-j + 2 .22 A 

+ 0.279 A X^-{ - 4.71 A f2-2.24APg2 - 0.284 Al^^^g 

. . (6,18) 

With a step load disturbance of 0.01 p.u. of pov/er 
in itrea 1, the responses of tlie state rariables I’e present iiog 
A f^ and A f 2 are detcnairEd and given in Figure 6.2. 
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6.3 TWO-ABSA U'C SYSTEM WITH BZCIT^TTIOil COKTSOI 

The block diagram of the tii70-area LBC system v/itii 
excitation control included in one of tiie areas is giren 
in Eigure 6,3. The voltage input for the second area, 
however, remairm the same as in Pigurc 6,1. Thus the 
effect of excitation control in one of the areas is 
studied in this section. 

• The salient differences between the block diagrams 

given in Eigure 6.1 and Eigure 6,3 are; 

(a) Excitation control loop included in Area 1 

The block diagi^m for the perturbation model of 

excitation loop has been derived from the standard IEEE 

30 51 

block diagram for type 1 excitation system ' • The 

same is shown in Eigure 6,4 in its complete shape. The 
description of tlie variables and parameters in t li i'w O c4X*.liO 
is 

A Ygi = texminal voltage 

Tg = regulator input transducer time constant 
K^,T^= amplifier gain and time constant 

Kg = constant corresponding to setting of the shunt 
field rheostat 
T-ti = exciter time constant 

Jzi 

K-i,T-,= gain and time coristant of the stafcilizei’. 

All the variables in the block diagram arc I'^erturbations 
from their nominal values. In this study the following 



1*sT 



Ft g i -3 Two-area system with exctt at 














simplifying assumptions are made: 

(i) = 0 

(ii) = 0 

(iii) the .limiter after the amplifier block is 
considered to be absent. 

(iv) the saturation function Sg is talcen to be a 

constant multiplier. However, the ramifi- 

cations of treating the saturation function Sg 
in its true shape arc dealt unth in Section 6.5. 

Also, to fit the above model in that given in 
I'igure 6.3 the following equalities are made: 

'"P 

1 

of S’igure 6,3 = |^ , of Figure 6,3 = ^ > 

'A |Yj| of Figure 6.3 

where Kj, Ig and ^ are described above. 


(t) Interaction of voltage on load demand 


In practical power sj^tems the effect of voltage 
on load forms a major interoctioa . Hence tMs is te-ken 
into account by means of the constants 8 


Y, 


and 


‘^21 represent the change of load per unit 

change in voltage in the Areas 1 and 2 resixjctively . A 


typical value of 1 ,0 is taken for each of these constants 
from the results obtained by Bogucki and hojcik'' ’.vho iiave 
conducted exhaustive tests on the Polish network for 


various kinds of loads 
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6.4 APii^IICATIOII OP OPTIMAL COI'JfROI THEORY 

The dynamics of the system shown in Figure 6,5 
can be given as 

x = Az+Bu+rAP^ (6.19) 

which, after defining the states vri.tii respect tc their 
steady state values, can be v/ritten as 

X = 1 X + B u (6.20) 

The state vector x in (6.20) consists of the 11 -elements 
/ Af^dt, A f^ 5 A , A Xg^-J ,/ A f2dt , A fg, A -t'g2 J * ^gv2 ^ 

A I Vi I , A V^^ and ^ last three variables are 

perturbed value of tcrmiixil voltage of Area 1 , exciter 
output voltage and stabilizer output respectively. The 


control vector u consists of the four-elements a ^c2’ 

^d A ,V 2 • Also 



.. ( 6 . 21 ) 
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Hiie nuaerical values for the exciter control 
parameters given in the A and 3 matrices are; 




1 .25 


= 0.05 ’./here 


12 = 0,05 r;here 


3p 


d1 


3P 


d2 


V 


3 r. 


^1 

* 


tiel 


a j - 

’hle1 


SjYgj 


3 ! ~ case may be. 


Pie Id tine constant 1^^^ 6,17 


Pie Id gain 


hid = 0-” 

0 *5 

Exciter time constant 1 ^ ^ - "10.0 

Exciter gain Kj 

Stabilizer time 

constant Ty = 1 ,0 




1 


-0.05 


= - 20.0 


Stabilizer gain ^ 0.04 

Aiapl:lf'i:-r gain Kj^ = 25.0; Amplifier time Oonstant 1_^=0.0 
Ihe rest of tlio parameters are as described in Apjendix A. 


The performance iiviex for minimization is taken as 

CO rp rr»j 

J = / (S H n + H u) dt (S.23) 

o - - 

Here Q. is a 12x12 matrix 'vVith tlie elements Q(2,2) and 
Q (6,6) each equal to 1, the rest of the elements being 
being zex*o. This means that veightages are given to the 
state variables A f-| , and Af 2 only. Also in this study, 
four diffei'ent values are used for E as detailed in 
Section 6.6. 


4 Por convenience in cciTiputstion tris gain is ine 
into the exciter gain making Xj=-500.0. Tliis 
by sliiftlng the line coming from S-, to the sum 
ahead of the amplifier block. 


orporated 
is dO'Tie 

a.lrg point 



Pontryagin’ s minimum principle is applied to find 
the closed loop controller is of the form 

u(t) = -K z(t) = -E"”’ 3^ P x(t) (6.24) 

where P (11x11) is the solution of the matrix Riccati 
e cfuation 

m *1 IT! 

(6.25) 
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P = -PA - A^P + P B R“'’ P - Q 


Phis matrix Eiccati equation is solved by the 
1112 

Blackburn -Hen ’ method given in Appendix C; however, the 

procedure given therein for determination of the initial 

solution talces considerable tine on the computer for systems 

of order 6 and above; hence the same cannot be used for 

the 11th order system studied here. Instead, the method of 
37 

Kleinman is used for finding the initial solution of tne 
matrix Eiccati equation. Piiis method in turn requires the 
computation of e"^"^ for which Yidyasagar ' algoritm is 

used. Both these algorithms are described in Appendix D. 
Phis procedure for detoimining the initial solution requires 
much less computer time as compelled to that given at the 
end of An OC- 1 V « 


6.5 SitPUHAPIOP I'DELIlSAftiPY 

Referring to Pigure 6.3, the block representing the 
satui'ation function Sg is a nonlinear function of As 

we are dealing with linear analysis, inclusion of the 
saturation nonlinearity in its actual form is not possible , 
Hoy/ever, as it is only a mild type of nonline arit3r, it can 
handled by the method of Gain Scheduling described belov;. 
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Gain Scheduling ; 

She methoi for calculating the saturation function 
Sg at a particular value of is shown in Figure 6.5A. 
Also let the saturation curve shown in Fig.6.5B be divided 
into, say, three equivalent linear portions OF, PQ and QR. 
It is easy to see that for each of the linear poi*tions the 
value of Sg ranains constant. Assuming a nominal value for 
Y^^, as long as the value for (7^^ + A ^ parti- 

cular linear range (Figure 6,5B), the corresponding value 
of Sg is used for computing the elements of the matrix A; 
and the matrix Riccati equation (6.25) is solved, Ihus the 
method reduces to that of solving three different matrix 
Riccati equations for the three values of S™ and these 

III 

tla’ee solutions, say Pg, P^ are used to determine the 
corresponding feedback gains IC-j , 1^2 and . YiAle deter- 
mining the response of the system prevision is made in the 
program to check at every step, in which of the linear 
portions the value of (Y^^ +A Y^^) is, and then the 
appropriate gain is used to compute the closed loop matrix 
(A - 3 K^)- i = 1,2,3. 





6.6 COMPUTATIONAL EESULTS 
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Four types of study are made in this chapter by 
employing four different combinations for the values of 
the parameter (= v/eighting 

matrix R, keeping the rest of the parameters constant 
in all the studies. 


a) 


b) 


c) 


d) 


3P. 


3P. 


3P, 




aP. 


ap- 




3 tv. 


3P. 


31-?. 


1^2 


3F. 


aiY. 


3F. 

a 


3! Y. 


3iY. 


=1.0, R = 


0.0, R = 


= 1 .0, R 


0.0, R 


1 .0 

O 

• 

O 

0.0 

~T 

0.0 

0.0 

1 .0 

c.o 

0.0 

0.0 

0.0 

0.1 

0.0 

0.0 

0.0 

0.0 

0.1 

1 ,0 

0.0 

0.0 

o .! 

0.0 

1 .0 

0.0 

0.0 

0.0 

0.0 

0.1 

0.0 

J 3,0 

o 

* 

O 

0.0 

0.1 

T.o 

0.0 

0.0 

0 . 0 ~ 

0.0 

1 .0 

0.0 

0.0 

0.0 

0.0 

0.01 

0.0 

^.0 

O 

• 

O 

0.0 

O.Oi, 

1.0 

0.0 

0.0 

0.0 

0.0 

1 .0 

0.0 

0.0 

o 

o 

O 

• 

O 

0.01 

c.o 

lo.o 

0.0 

0.0 

0.01 


The solution matrices (symmetric) P of the 
matrix Riccati equation (6,25) in the four cases 
respectively are 
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1352 - 0.1078 - 0.1219 - 0.0254 - 0.1352 0.0044 0.0078 0.0017 0.1867 0.0017 - 0.0040 
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0,3639 - 0.0236 - 0.0947 - 0.0260 - 0.3639 0.0236 0.0947 0.0260 0.0295 0.0010 - 0.0109 
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In each of the four cases above, the closed loop 
—1 5 } 

system matrix (A-B R B I) is computed. The responses 
of the frequency deviations Af-j and a £2 obtained for the 
cases (a) and (b) are given in Figure 6.6 as continuous 
and dashed curves lespectively . Also the responses obtained 
for the cases (c) and (d) are given in Figure 6.7 as 
continuous and dashed curves respectively. 

6.7 COlCIUSIOhS 

The effect of excitation control in one area of a 
tv/o-area IFC system is studied. Figures 6,6 and 6.7 summarize 
the results obtained. The system studied by Burick and 
responses obtained therein are given in Figure 6,2 for a 
comparative study. 

The follov/ing conclusions are drawn from the 
responses obtained. The interaction of voltage on load 
demand prevents the system from reaching at least a non- 
zero steady state and thus has a deteriorating effect on the 
system response, Becroasing the penalty factors on the 
two inputs cor*responding to exciter control in Area 1 and 
voltage perturbation in Area 2 will prove to be much 
harmful to the system. 

In the study made with the assumption that there 
is no interaction of voltage on load demand, the steady 
state frequency error is lesser than that obtained in 
Burick’ s study (Figure 6,2) This is explainable to be 



time, sec 






time,sec 
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due to inclusion in this study, of tiie usual load 

frequency controllers AP - and Ap „ 

Cl ^ c2 

A comparison of the responses given in 
Pigux'‘es 6.6 and 6,7 with those obtained by Elgerd and 

3 

Posha emphasize the optimistic nature of the results 
obtained with the assumption of noninteraction betr/een 
Megawatt frequency and Megawar voltage control loops. 
The above results also confirm that in IMJ systems 
where the inherent stability of the system (without 
considering exciter loop) is only maiginal, inclusion 
of exciter loop will actually lead to instability of 
the system. Thus the study made in this chapter, is 
great importance because it indicates the necessity 
or othemvise, of employing compensating networks that 
will stabilize the system. 



OHAFPER 7 


ROmiEASURiLBILITY OF STATES AIID COIISTSUCTIOF 
OP AIJ 02SEEVER 

Tbs optimal closed loop controller designed by 

3 

Elgerd and Posha for a two area system is based on the 
assumption tiiat all the state variables are available for 
feedback; In a practical system measurement of some of 
the state variables may be either extremely difficult or 
altogether impossible because of limitations in tlie measur- 
ing apparatus. One method of overcoming such a difficulty 
is, constructing an observer of the luenberger type ’ 

By this means, an observer is constructed for the in- 
accessible state variables which are then used in tlie 

optimum regulator configuration for finding the responses 

8 

of the state variables. Bongiorno and Youla have given 
the theory of construction of a compatible obser’/'er of 
the Luenberger type. Section 7.2 gives a summary of 
results obtained by them. Section 7.3 deals with the 
computational results for the construction of an observer 
for the two-area LPC system. 

7,2 THEORY OP COHSTHUCSIOH OP A OQYPATIBIE OBSERVER 
(LUENBERGER TYPE) 

Let the system dyoamic equation and output equation 
be given by 

X = A X + B u 


(7.1) 
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and 

Z = D z (7.2) 

Here x is an n- vector 
is a p-vector 
u is an r-vector. 

A is an (nzn) system matrix: 

B is an (nxr) input distribution laatrix 
D is a (pxn) measurement matrix;, 

D is constructed to have a simple structure; each 
row consists of 1 in only one place corresponding to the 
observable state and the rest of the elements in the ro’v 
are zero . 

Por this system, the dynamic observer is defined 
as 

y;=Hv + Sz + Eu (7.3) 

v/here v is an (n-p) -vector 

H is an ( n-p )x( n-p ) -matrix 
S is an (n-p) x p matrix 
E is an (n-p) x r matrix, 

H is chosen to be a stable matri:i with eigen- 
values different from those of the system matrix A of 

(7.1). 

Using this observer, the estimate x of the 
state vector x is given as 
X = L-| z + 1(2 Z 


(7.4) 
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where is ah nxp matrix 

i 

112 is an iix(n-p) matrix. 

Definition ; 

An observer is said to be compatible if its output 
equals the state of the plant, to within an oxponcntially 
decaying error. 


For the observer to be compatible, it is necessary 
that (a) the plant be completely observable and (b) the 
observer be completely controllable with respect to tlie 

plant output. Matnematically, condition (a) means that the 
matrix 

Di . .1 . .,2.2 2 m n-l rjT) 

, .AD U ) D : : (A ) (7.5) 

should be of rank nj and condition (b) means t.iat the 
matrix 

= s : H s : s : . . 

• • « 

should be of rank (n-p). 


A = 
o 


. tP^-I 

. h 3 
_J 


(7.6) 


2he observation vector y is related to the state 
vector X by 


(7.7) 


y = U X + e 

where U is tne unique solution of 

D A - E IT = S B 

provided A and n do not have common eigenvalues. For a 
choice of E = U B, tlie vector e in (7.7) is an exponentially 
decaying error between the actual and observed staae vectors, 



(7.9) 


The differential equation for tliia error is given bv 
£(t) = H e(t) 

The solution of the equation (7.9) is 
e(t) = exp (H(t-tQ))e(tQ) 

The reconstructed state vector is nov/ given by 
X = I(f X + £) 
where 


of l 

1 



L^J 


exp(H(t-t^))e(t^) 


(7.10) 


(7.11) 


(7.12) 


0 ^ in the above equation denotes a p-element zero vector. 
The matrices 1 and W can be split up as 


I = |li :iq 

and 



U = ^ (H) U 

and 


(7.13) 

(7.H) 

(7.15) 


if = Q A 

c 0 

where and r 

V 

tively. Also 

f 

1=0 

with is 

A and ’s are 

X 


(7.16) 

^ 8j?e de±iix8G in (7*5) and (7*6) rsspec— 

0 

1 

“i 1 (7.17) 

the characteristic equition of the matri^c 
its coefficients. The siatriii Q in (7.16) 


is defined as 
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“1 4 

“2lp 

“ T 

n-1 

a T 

n ■%) 


“2 Ip 

»3 ip ... 

ct I 

n -p 

-pp 


♦ 

« 



« 

• 

(7.18) 

°n-l 

“ n -p • • * 

%P 

-pp 


“n ^ 

Qpp ... 

0 

-pP 

0 

-pp 



Here ^ and denote pxp identity natrix ard pxp null 
matrix respectively. 

For the observer to be compatible, it is required 
that the matrix W should be nonsingular. With this condi- 
tion, 1 is computed in the simplest manner as equal to 
and the estimate x of the vector z computed as in (7.4). 

7.3 APPLIClflOIv TO A LOAD PHSQUE'JCT GOIIiHOl SYSIEI" - 

CClIPIJTATiaiAL ESSUSDS 

Figure 7.1 gives the block diagraa of a t\’7o-area 
iFC system. The description of the matrices A (9x9) and 
B (9x2) of (7.1) is as given in (A-8) ana (A-9) of 
Appendi^r Aj also they are computed as given Section 3.7 
of Chapter 3. Ene matrix D (4x9) of (7.2) is taken as 

1 0 0 0 0 0 0 0 0 ~ 

0 1 0 0 0 0 0 0 0 

0001 000 0 0 
OOOOOOlOO 


B = 
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Thus only tho four state variables and x^ are 

considered to bo measurable . An observer is to be cons- 
tructed for this system from the four available outputs 
and the two inputs. Hence the dynamics of the obseiver 
should be of fifth order. As stated above, the matrices 
H and S of the observer are so chosen that they form a 
controllable pair(H,S) and also that the matrix H does 
not have any eigenvalue in common with the original sys- 
tem matrix A, Keeping these considerations in view the 


matrices H and 

S are 

chosen 

as 







■^1 ,72 

0 

0 

0 

0 ~ 


T 


0 

0 


0 

-1 .87 

0 

0 

0 


0 

1 

0 

0 

H = 

0 

0 

-1 .99 

0 

0 

S= 

0 

0 

1 

0 


0 

0 

0 

2.15 

0 


0 

0 

0 

1 


0 

0 

0 

0 

-2.33 


1 

0 

0 

0 


The coefficients of the characteristic equation (7.17) 
are computed and the matrix U of (7.16) is also computed, 
Erom (7.15) u is computed as 


0.581 

-1.539 

-0.356 

-1 .325 

o 

* 

o 

1 

1.539 0.356 1.325 

0.410 

0.0 

1 .213 

0.249 

1 .020 

0.320 

-0.678 -0.139 -0.570 

-0.179 

0,0 

0.939 

0.399 

1 .037 

0.329 

-0.939 -0.173 -0.771 

-0.245 

0.0 

0.452 

0.071 

0.358 

0.115 

-0.452 -0.227 -1.150 

-0.370 

0.429 

-0.155 

-0.020 

-0.117 

-0.038 

0,155 0.020 0.117 

0.038 
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I'll© matrix E is computed from the relation E 


-5.121 

5.121 

3.997 

-2.234 

4.111 

-3.058 

1 .441 

-4.628 

-0.480 

0.480 


The matrix Vf is also calculated as in (7.14); and finally 

•r , -1 

1 Which is the inverse of 7 is computed. Eor use in 


up as L-j 

and Lg . 

Here 


1 .0 

0.0 

o 

* 

o 

0.0 

0.0 

1.0 

0.0 

O 

• 

O 

■6.044 

2.046 

0.020 

-0.005 

0.0 

0.0 

1 .0 

0.0 


Li = 


0.008 -5.330 -3.204 -0.001 

■ 3.966 - 0.022 - 0.022 - 0.002 

0,0 0.0 0.0 1.0 
55.160 2.403 1.005 -1 .977 


1 

L 

-105.479 

-7.482 

-3.121 

5.527 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.815 

-3 .825 

6.267 

-0.045 

12.979 

0.0 

0.0 

0.0 

0.0 

0.0 

1 .668 

9.968 

-4.836 

-0.008 

-2.278 

0.809 

1 .912 

-0.242 

0.977 

8.144 

0.0 

C.O 

0.0 

0.0 

0.0 

16.217 

-4.493 

11 .925 

1 .455 

-103.892 

■50.668 

15.991 

-37.043 

-8.423 

314.404 
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The dynamical equations of the original system (7.1) 
and the observer system (7.3) are then solved for obtaining 

A 

the response v/ith inputs u-j and U2 given in tcims of x of 
(7.4). An initial condition x^(o) = -O.OI is taken for the 
original system and a corresponding one for the observer. 

The coefficients for the optimal controller structure are 
computed by applying the Pontiyagin’s minimum principle 
(Appendix B) to the system under consideration. As the 
system taken in this study is identical to that tairen by 
Elgerd and Posha^, the feedback coefficients computed by 
them for the optimal controller structure, are utilized 
here . 

In the two -are a IPO system the state variables of 
interest are the frequency deviations in the two areas A f -j 
and Afg and the tieline interchange error A 

responses of these state variables are determined (Pigure 7.2) 
by cascadirjtg the observer system with the original system. 

Pigure 7.3 gives the response of the above state variables 

3 

under the assumption of complete state feedback . 

7.4 COlvESlTS Ai'ID COITCIUSIOIIS 

In this chapter an observer of the Inenbcrgcr type 
is designed for tlie tv/o-area IPG system assuming that only 
four out of the nine state variables are measurable. The 
estimate of tlie state vector thus obtained is utilized for 
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determining the response of the system; the response with 
complete meas’orement is also giTen for comparison. 

It is seen from the Figures 7.2 and 7.3 that the 
responses obtained by using an observer are almost iden- 
tical to those obtained by assioming complete state feed- 
back, and hence tiie method is advantageous in load 
frequency control of multi-area systems. However, when 
the observer system is of high order its cost will be 
proportionally high. Under such a situation, it is 

better to go in for some kind of suboptimal control. using 

4l 

the available state variables for feedback purposes 

The theory of observers ensures the reconstruction . 
of the state vector v/ithin a small error under certain 
conditions. It is desirable to keep the Gigenvalues of 
the observer system dynamics, to be reasonably large so 
as to er^ure fast decay of any possible error in the 
initial condition telnn. A second requirement is ehat the 
observer should have smaller eigenvalues in order that its 
hardware realization is feasible. A suitable compromise 
has to be made be fcvoen these t7/o conflicting requirements, 
Though the theory for the construction of an observer 
is fairly complete and is applicable to IilO systems, a more 
rigorous treatment, which also takes into account, the 



physical realities, will be, to consider noiss in the 
system djaiamics as well as in measurement, Por optimally 
controlling suen linear stochastic IPG systems, well 
established s i^ocnastic control techniques are available^^ 
and have sinco been applied to small order IPO 3,'/stem^^. 



CH/iEPSE 8 


C0IICIUSI0IT3 Aim 3U&G3STI0irS 

State variable approach and optisal control tech- 
niques aro nov7 extensively used for the analysis of ITO 
systems. Hov/over, when they arc applied to higher order 
systems, computational effort and time increase considera- 
bly. The methods pixisented in this thesis, for the con- 
struction of suboptinal controllers, 7d.ll enable one to 
overcome this ■ difficulty and thus v/ill enhance the use- 
fulness of optimal control tochniques in the area of Load 
Ere quo ncy C ontrol . 

The methods given for the model reduction of linear 
IPO systems arc both elegant and efficient. In the ease 
of suboptimal conti'ol by aggregation, it is also possible 
to correctly and quickly check the accuracy of representa- 
tion by computing the degradation in pcrfomanco, 

25 

Kiniesy and Bohn- have constructed a suboptiitial 
controlloi’ having a tv/o level structure for the large signal 
model of a multiaroa liPC system. Por the first time, an 
attempt is mado here for the reduction of the las'go signal 
model of a multiarc a LPO system to a lov;er order nonlinear 
model and also for the construction of a suboptimal control- 
ler in the shajX; of the optimal controller, dcts.rmined by 
Lukes' method, for tlio reduced model* 
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The studies made here by including excitation 
control loop in IPO systems not only reveal the behaviour 
of the system uiidei- realistic conditions but they also 
confirm the liiglily optimistic nature of the results 
obtained by the assimiption of noninteraction of tiio 


Megawatt frequency and Megavar voltage control loops . 

As such tlTcso studies arc v^ry helpful in tAiing appro- 
priate decisions regarding provision of cosipcnsating 
nctv/orks for systcras in which the open loop system matrix, 
without the inclusion of exciter loop, is marginally 


stable , 


Tlie I'csults obtained with the construction of an 
observer for an IPC system arc quite satisf actoiy. However, 
the cost of construction of an observer increases in pro- 
portion with the order of the observer system which in turn 
depends on the order of the laultiarca system^ as such the 
utility of the method is confined to small order multiarea 
systems only. 


Some of the related problems for further re 


ere] 


aro: 

(i) Oi^tiEial deterministic control of multi-area LPC 
systems with thn inclusion of both tlx governor and 
ticlino power nonlinear i ties. 

(ii) Seeking reduction techniques for tlx optimal 
control of LlfS sj^stems, consisting of many rnxas (say 
6 and above ) which arc different in their sis-o (if/ 
capacity) and/or paramctvrs. 
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(iii) Extending the theory of aggreg-tion du;. to Aoki 
for the rv^duction in computation of cptimel controllero 
for linoi.r soocnc.otic uiultiorca EEC systerns . 

(iv) bubo2)tiraal stochastic control of the lar^^o 

O 

signal modol of an IiPC system having ticlinc power 
n onl inc arity . 
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appeioix 


IvrO-nuBEA iOin) PBBQUEMJY COIKIROI SYSTEM 


T} 3 c ty/o-arca load fa^cquoiacy control system considered 
by llgcrd and Poaha^ is taken up for study. The block 
diagrcaa for the two-aroa system is given in Figure 7 . 1 . The 
system differential equations arc; 


dt^-^ *^tie1 ~ ^ 12 (^ Af^dt - /fif^dt) 

|t( /A%dt) = ^ f^ 

dl "^i + I^i ^ fi + A P^^-a 




1 


at ^gvi= - T^-'^gvi - 0^^‘fi t H 


1 


■gvi 


&p 


ci 


&p 


tio 2 


^ 2 * ^tiol 


p"- ^- = ’>2 

^ ‘ ^tiol 


(■ 1»1 ) 

U.2) 

U. 3 ) 
(A , 4 ) 

(ii. 5 ) 


In viow of this relation f a P »^+ 4.^ j. 

a.xion, / & P^j_Q2 3 .S not oonsidered 

as a state variable. 


^2 “ '' 2 > 


«* - s* 

12 ~ 

f* 


300 


6 0 Hz 

4 . 

^1 ~ ^ ^ ^ ' 


“ ®gv 2 “ 

= Ig = 2.4 Hz/p.u, m;7 

■# 

T ^2 w 0.545 p.u. Mi? 

'^r 1 ** ^r 2 * 
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- 2)2 = 8.33x10*"^ p.u.Lf?/H 2 A = 0.01 p.u. IIP 

’^tl = ®t2 = 0-3 soc. s ^ 0.0 


Equations (n.1) to {A. 5) can be put in the fern of vector 

dif icroiitial equation as 

X = A X + B u + fa P r ev 

where x is a 9 olcmont rector eensisting of the state 
variables / aP^ . dt , / Af^dt, a f , a P^^ , A , / A f^dt , 

Af^j '^^g 2 ^^grP’ ]i is a P-elcnent rector consisting 

of the control rariablcs A p^^ and A p^^- ^ 3 o a p^ is a 

2-oloiii9nt rector consisting of the two load disturbances 
AP^^ in Airoa 1 and A 2 as its elcnents. r is 

ci load disturbance distribution matrix given as 
f*. 


r = 


0 0 - 


2H 


i 


0 0 


0 


0 0 0 


0 C 0 


f* 

OTT 

eng 


0 0 


0 0 


By rodefiniiig tlie state and control rectors in toms 
of their steady state values, the terms r A P^ in (,..6) ecai 
be eliminated and the equation (ul.e) c- b- ... 


2 X “f- 


. 7 ) 
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Horo A = 

0 
0 


0 


0 


0 

0 

0 

0 


■12 


0 




12 


0 

• 1 


2Hh 


0 


0 


0 

-1 


0 
0 

1 


1 2H^ 


=1 

T 


t1 

0 


0 


0 

0 

0 


-I 

^gvT ' 


-^ 12^^12 


0 


0 0 0 


0 0 


0 0 


0 ■ 0 


0 


0 


_fn 


12 


f^T 


0 

* 

12 


2Hi 


0 


0 


0 


0 0 
0 0 

0 0 

0 0 

0 0 


!12£!32 -£!^ ^ 

2 B 2 ' - 2E^ ■ 2^2 


0 


0 


0 


-1 


0 

0 

0 


0 

0 

0 

1 




n2 ^t2 
-1 ^ -1 




gV2 


T 


gv2 


. . (A.8) 


B= 


0 0 0 0 


T 


gvi 


0 0 0 


00000 00 c 


0 

1 


r 


gv2 


T 


(A.9) 


Por applying optiiaal control theory to the above system 
the perfoxmaneo to be minimized is tcAcen as 


CD 


a? 


<J = / (x Q s. a ^ H) 


(... 10 ) 



iiP£EItt)IX B 


(mimjj COiSEOLIiBR USING P03:5TRY.^GIN' S 


HINIMuli 


priiiciple'^'^ 


Let tile system dynamics be of the fora 


i; = A X + B u 


0.1) 


Here x is an n-vcctor, u is an r-vector, A is an nxn matrix 
and B is an nxr matrix. 

The cost to be minimized is toJeen as 

Cq 

J = ^ / (x^ Q X + up E u) dt (3.2) 

0 

where Q, and E arc nxn and rxr matrices respectively. T/c 

‘T 

augment the cost J by adding x"* (-‘'•x + Bu - x) v/hcrc ^ is 
an nxl costatc vector. 


J' 



1 m jj 

^(x^ Q X + u R u) + X. 2. + B u 


±)j d t 

4b.3) 


Let the Hamiltonian be defined as 

H = ^(x^ Q 2 + B u) + X + B u) (3.4) 

Also let the costatc vector \ satisfy the differential 


e quation 

1 = - ^( 2 , L, 2 ) Cb.5) 

Equation (B,1) can also bo written as 

3H 
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Lot U "be a set of all admissible control vectors 
and u* one such admissible control, the correspon- 

ding solution of (3.5). Then Pontryagin's minimum principle 
says that for u* to minimize H, 


= H(x*, u*, X*) = min H(z*, u, (B.7) 

UeU 

y/herc is the scluticn of (3,1) with u--' substituted 
the min, for u. Also minimization of H in the above 
manner is required for the mirnmization of J . A necessary 
condition for minimizing H as given by (3,4) is 


-J.K. = 0 (B.8) 

SU' , 

On performing the differentiation given by (B.8), we get, 


I 

R u + B X = 0 


or 


* 

u_ 




(3.9) 

(3.10) 


Using the value of u* given in (B,10), (B,5) and (3.6) 
give the differential c quationsthat 3^ and must satisfy 

X* = A X* - 3 B^ (B.11) 

- Q X* (3,12) 

V/c v/ould like u <s.s a function of x| so wo assume 

X = P X (B.I3) 


and try to determine the conditions on P. Substitution 
for the costate in terms of x into (B.II) and (B,12) gives 
p 2;* + P X* = -A^ B X* — (i X* (xi,14) 
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X* z= A X* - B R"*"^ B X* (B.I 5 ) 

Substituting x* into (B.H) and factoring out on the right 
the rector x*, gives the nxn matrix differential equation 

P = -P A - P + P B R“'’ B^ P - Q (B.16) 

This matrix differential equation is called the 
matric Riccati equation. Por tlx infinite tiao problem, 
the Riccati equation has a steady state solution f^r which 

p = 0 (B.17) 

Once the steady state solution to the matrix Riccati 

equation is found 

-1 T -1 T 

u* = -R B L = -R B P X (3.18) 

and if we let 

K = R-’^ B^ P (B.I 9 ) 

then 

u*= -K X (B,20) 

which is in the desired ferro. 



iiPPEiroix 0 


SOIUnOlT OE IE-J11E.I1 RlCCaill BQU^J?I01'T BY liBV/TOIJ-IL'^KSOil IffilHOD 

(BLAOKBUSl^ - MM ‘iliGORIK-M) 


BiiC matrix differential equation that one coacs across 
in the application of optimal control theory to dynamical sys- 
tems is given by 

P = -PA - ? + P 3 P - <3 (C.1) 

Hero A is the system open loop matrix, B is the control 
distribution matrix, Q weights the state variables and E 
weights the control variables. Both Q and R occur in the 
performance integral. P is the matrix fer 7/hich solution is 
sought for. In the ease of an infinite time problem, (l.T) 
reduces to the natrmc quadratic equation given by 

0 = -P A - P + P B R”"' B^ P - Q {C.2) 


Many methods exist in mod' m control 


literature for 


the solution of this 
1 1 

them Blackburn has 


matrix quadratic equation. In one of 
applied the Rewten-Raphson method for 


the iterative solutiun of (C,2). Ills iterative sanomo is 


essentially as follows: 



(. - sp^)Ai 

(0.3) 


where G = 

iteration, 

tiVGly,a 


1 I 

B R*" B and the subscript k indicates the kth 


and the notations 




and { D } denote, respcc- 


X n(n+1) dimensional co3.umn-vcctor and 



ITS 


iiiirolificd the 


1 1 

2 n(n+ 1 )x ^n(n+l) raatrix formulated from the rixn matrix 0 

~ 2 2 12 
and the n xn matrix D, Man has further 

above method as follows: 

Writing (C.3) as 


■k+1 


= [h] 


+ { Ix ( ) ^ + ( ) ^xl ^ 


-1 


(0 .4) 


Begrouping th.c terms in (Cv4) give 

-1 


pk+^ = +(^--GPk) XI > + [PfcJ 

- {lx(A-GPj^)®+(h-GPj^)^xl} ~ |Pg-(A-^?j^)+(;-GP3^)^?~- 

•* (0.5) 

V/ith the use of the Kroncckcr product, it can be shov/n 
that 


\(A-GPk)+(^.-GPj.)^ = {Ix(a-GP3^)^+U-GPj,)IxI1 


.. ( 0 . 6 ) 


Substituting (0.6) in (0.5) yields 


-1 


"k+l 


= { Ix(A-GPj^)^ + (^-GPj^) xl> [-Pi;(JPk-<il (0-7) 
In view of (0,6), (0.7) can be further simplified 


to give the desired iterative scheme 


( 0 . 8 ) 


Equation (0.8) is nothing but tto matrix lyapuncv aquation 


of the type 

S P + P^ S = -P 


(0.9) 


for which many methods exist 
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One problem still remains viz. finding the initial 
solution for the it:z-ativc scheme given by (0.8). 'The 
procedure usually followed is to solve the ^n(n+1 ) simul- 
taneous scalar differential equations resulting out of 

the raatrix differential equations given by (C.l) till such 

20 ' 

time that the resulting matrix P stabilizes the closed-loop 
matrix (h-B E B P) . The matrix P thus got is used as 
the initial solution for the equation given by (0,8). 



APiESDIX D 


PIITDBTG lUITIBi SOLUTION FOR Ilr-THIX RICCLTI EQBBION 


It is well known thot tix coiiTcrgencc of tho 
iterative process given by (0.8) of Appendix 0 is assured 

if the initial solution matrix stabilizes tlx closed 

—1 2 

loop system matrix . (A-BE B P^). 

37 

Kleinman has shovrn that if the original sysrtem 
given by 

x=Ax+Bu (B.l) 

is controllable, then the raatrix vV”"' v/ill stabilize tho 
matrix A whore 

Y{ - ^ B B^ 0 ”’““’ dT , T = arbitraiy (B*2) 

o 

• T-1 

Thus the appropriate initial solution will nov/ be Pq=-*" 

Further, in (D,2) tho matrix exponential e"‘ can be computed 

36 

in a closed form by the algorithm given by Vidyasagar 
which is given belov/. 


^ t 

Determination of e ~ in a Closed Form 


We can vorite e as 


At v’ a. .i 

e = I 1 xi 

i=o 


0 . 3 ) 


where i = o, n-l are functions of t. It is the 

aim of the algorithm given by Vidyasagar to determine these 



scalar functions. The algorithm runs as below: 


(i) find the cigcnvaluos cf the matrix 1, 

* n 


(ii) Compute (“^) n ( X. - X. ) 


3 


3 1 ' 


where P(a) = H (s- ^ • ) = I a- s' 

« Ji. • <ir 

1=1 1=0 


(D.4) 


(I).5) 


(iii)De fine 


n r e 


,(t) = I x^r { 


0, . . . ,n-1 (P.S) 


^ 3S^ X 


Then the scalar functions are dctc rained by 


(iv) aj^ 


(D.7) 


i — 1 5 * ♦ * • u 


(v) Pinally e'^'*' is computed from (P.3). 



ilPPEEDH E 


(i) M atricGS and Vectors Described in Ohaptor 4 

The elonients of the right hand rector rj_ (20:c1 ) cf 
( 4 . 69 ) arc described belov/. Eat 
■f-j 2 f*x0,5x / 

p ^ __ __ 

Then 


n (1 ) = -2pP-j 2 

T\ 

(11) = 0.0 

n ( 2 ) = - 2 P P22 

n 

(12) = 0.0 

n (3) = 0.0 

n 

(13) = 0.0 

( 4 ) = -'2pP23 

n 

(14) = 0.0 

n ( 5 ) = 0.0 

n 

(15) = 0.0 

r, (6) = -2PP24 

n 

( 16 ) = 0.0 

( 7 ) = 0.0 

n 

(17) = 0.0 

( 8 ) = 0.0 

T 1 

( 18 ) = 0.0 

11 

0 

• 

0 

r\ 

(19) = 0.0 

(io)= 0.0 

n 

( 20 ) = 0.0 


Here P^-^'s are the elements of the matrix P* commuted as 
13 * - 

under ( 4 . 64 ). 


How lot = 2a^ + • 

^ijk“ ^ii 


i = 1,2, 3, 4 
3=1 ,2,3,4 
k = 1,2, 5, 4 


Also let the elements a. - described below to be the 

-L J 

elements of the mteix A* = (^ + BK) computed as under 
( 4 . 66 ). Then the matrix (20x20) of (4.69) is doscribed 


as 
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Tho elements of tbe right liand vector e (35x1) 
aro computed as follows: Let 

0,5f* '^12 ©.SSSf*" 

t = f7 and t_ 

let 


11 

J’ 


C 12 

- 

'=16) 



If - 

ty(c 

32 -c 

35 ) 

% = 


°12 

- 

*=16) 



% 

t^(c 

42 " 

°46) 

0 

11 


°22 

- 




% = 

ty{C 

42 ■“ 


= 


^22 

- 

<=26) 



If = 

t^(o 

52 " 


<30 = 


C 32 


^ 35 ) 




ty (0 

52 “ 

^=56) 

Here 

"id 

' s a 

iTG 

; the 

e 

lemonts 

of the ■ 

transfoimation 

( 5 x 8 ) 

' . 

Also 

let 







=’ 

P 11 

^a 


Pi 2 

<30 

+ PI 3 

■5c + PI 4 

% 

P 15 

% 

II 

P 12 


+ 

P 22 

'3c 

P 23 

■^e + P24 

■^g 

^25 

If 


fc - P13 'la P 23 ^ ^^33 ^°34 % %5 % 

% = P14 + P24 % + P34 + P44 % + P45 % 

fe = P15 + P25 'ic + P35 % + P45 % + P55 


whore p..*s aro the Glements of the 
matrix Rice at i equation. Also lot 



in^'s be its elements. Let 


?/-[ 

= 3111 

*" 3^51 

I 

rr 

II 

-%4 

W2 

= hi2 

- 3I52 

W5 = hi5 - 

■ ^55 

Wj 

= 11^3 

- 3153 
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Then ^ (1 ) = 

e(2) = 4w-|W2fg_ + 27/^f^ 

^(3) = 2v4f^ + 
e (4) = 4w-,w^fg_ + 2wff^ 

^ (5) = 2Vir^f^ + 4v7-j^'^^xc 

e (6) ^ 4w^w^f^ + 2wff^ 
e(7) = 2w| + 4w^w^f^ 

e (8) = 4w-jV/^f^+2w^f^ 
^(9) = 2w^ 
e(lO) = 2v|f^ 

E ( 1 1 ) =: 4w27'^jX'-j^+2v;2fQ 

e(12) = 2w|f^ + 


z ( 13 ) = 4W274^f^+27/^ 

e (14) = 2vT^f^+4W2W^f^ 

A Q 

E (15) = 4w2W^f^+2w2fQ 
E (16) = 2w^f^+4w27/^f^ 

E (17) = 2w3f^ 

^ (18) = 4w3W_^f^+27|f^ 

£ ( 19 ) = 2vy| + 4w3V7^f^ 
e (20) = 47/37/3 f^+2w3f^ 
e (21) = 2r/3f^ +47/3^3 
e (22) = 2w? 

e (23) = 4w^W3f^+ 2wjf^ 
e (24) = 2w3f^ + 4w^?/3f^ 


e (25) = 2w|f^ 

e (26) = 4v72V/3f^+4w^W3f^+4w^V72f^ 

E (27) =. 4 v 72 'v 7 ^f^ + 4 w^w^f^+ 4 w^W 2 f^ 

e (28) = '*’ ^^"''^l''^2^e 

e (29) = 4 w 3 W^f^+ 4 w^w^ + 4 w^W 3 f^ 

E (30) = 4^3 W3f^ + 477^^3^^ + 4w^773f^ 

E (31) = 477^W3f^ + 4w^773f^ + 4T/^^^fo 

A 

e (32) = 4v/37/^f^ + 4w2W^f^ + 4772^/3^^ 
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^ (34) = ^''^''’^2^'^4^o 

e (35) = 4w,Wj-f + 4w^v/ f + 4-^,;v,f 

■-i- J ^ ^ D (1 J O 


Hhe elements of the right hai^id vector 


0 OI:ipU"tCu 


o. 

w40 


e(l0) = 0.2G673xlo'’ 
e( 1 l) = 0.33333x1 o”' 
e( 12) = -0.55140x10'* 

E(17) = -0.69565x10'* 

Iho rest of the elements being zero. Also the clcncnfe of i>li 

solution vector u ore computed os 
n(lO) = -0. -I- 43 60x10*^ 
ja(11) = 0. 56631x1 0“"^ 

m(l2j = -0.12688x1c'’ 

m(l7) = 0.96424x10° 


The rest of the elements being aero. 
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